Chiral three-nucleon force at N 4 LO II: Intermediate-range contributions 
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We derive the subleading contributions to the two-pion-one-pion exchange and ring three-nucleon 
force topologies emerging at next-to-next-to-next-to-next-to-leading order in chiral effective field 
theory. The resulting expressions do not involve any unknown parameters. To study convergence 
of the chiral expansion we work out the most general operator structure of a local isospin-invariant 
three-nucleon force. Using the resulting operator basis with 22 independent structures, we compare 
the strength of the corresponding potentials in configuration space for individual topologies at 
various orders in the chiral expansion. As expected, the subleading contributions from the two-pion- 
one-pion-exchange and ring diagrams are large which can be understood in terms of intermediate 
' excitation of the A(1232) isobar. 
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I. INTRODUCTION 

73 ' 

Three-nucleon forces (3NF) are presently subject to intense research, see Refs. (ll-fioj for a selection of recent few- and 
i " i . many-body calculations along these lines and Refs. [HI [13] for review articles. On the one hand, rapidly increasing 
computational resources coupled with sophisticated few- and many-body methods allow nowadays for reliable and 
t-H , accurate nuclear structure calculations for light and even medium-mass nuclei. One can, therefore, relate the properties 
^ • of the nuclear Hamiltonian to observables in a reliable way and without invoking any uncontrollable approximations. 
On the other hand, considerable progress has also been reached towards quantitative description of nuclear forces 
using the framework of chiral effective field theory (EFT) , see recent review articles [r3 - fl5j and references therein. In 
' particular, nucleon-nucleon (NN) potentials at next-to-next-to-next-to-leading order (N 3 LO) in the chiral expansion 
were developed [l6|, [l?} , which allow for an accurate description of NN scattering data up to laboratory energies of 
the order of £i a b ~ 200 McV. For heavier systems, the accuracy of theoretical predictions is currently limited by the 
\ 3NFs for which only the dominant contributions at ncxt-to-ncxt-to-leading order (N 2 LO) in the chiral expansion of 
I ■ the nuclear Hamiltonian have so far been employed in few- and many-body calculations. 



. — ■ The cn i ra l expansion of the 3NF at one- loop level, i.e. up to ncxt-to- next- to- next-to-next-to- leading order (N 4 LO), 
k> ' can be described in terms of six topologies depicted in Fig. [T] The first nonvanishing contributions to the 3NF emerge 
at N 2 LO from tree-level diagrams corres pon ding to the 27r-exchange, one-pion-exchange-contact and purely contact 
graphs (a), (d) and (f), respectively [Til. Il9j|. The shorter-range terms emerging from diagrams (d) and (f) depend 
on one unknown low-energy constant (LEC) each which can be determined from suitable fcw-nucleon observables, 
see e.g. [1, [fl EH, H(| ■ The long-range contribution (a) is, on the other hand, parameter-free since the LECs ci, C3 
and C4 accompanying the subleading ttt:NN vertices can be extracted from pion-nucleon scattering, see [21rl24| for 
heavy-baryon results, Refs. [25l.[26j for some more recent calculations using manifestly covariant formulations of chiral 
perturbation theory as well as Refs. [27j for an attempt to determine these LECs from proton-proton and neutron- 
proton partial wave analyses. The resulting 3NF at N 2 LO has been intensively explored in three- and four-nucleon 
scattering calculations, see [ll| and references therein. One finds a good description of low-energy nucleon-deuteron 
scattering observables except for the well-known, long-standing puzzles such as the vector analyzing power in clastic 
nucleon-deuteron scattering (the so-called A y -puzzle) and the cross section in the space-star breakup configuration, 
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FIG. 1: Various topologies contributing to the 3NF up-to-and-including N 4 LO: two-pion (2ir) exchange (a), two-pion-one-pion 
(2n-ln) exchange (b), ring (c), one-pion-exchange-contact (d), two-pion-exchange-contact (e) and purely contact (f) diagrams. 
Solid and dashed lines represent nucleons and pions, respectively. Shaded blobs represent the corresponding amplitudes. 



see Ref. for more details. Promising results for low-energy four-nucleon scattering observables based on the chiral 
3NF, especially in connection with the A^-puzzlc, are reported in Ref. Q. While nucleon-deuteron scattering data 
at higher energies are also well described, the theoretical uncertainty increases rapidly reflecting similar pattern in 
the two-nucleon sector at this order in the chiral expansion. Promising results based on chiral nuclear forces were 
also obtained in nuclear structure calculations showing, in particular, sensitivity to the individual terms of the N 2 LO 
3NF, see [ll| and references therein. All these findings clearly underline the need to include corrections to the 3NF 
beyond the leading terms at N 2 LO. 

The first corrections to the 3NF emerge at N 3 LO from all possible one-loop diagrams of type (a)-(c) involving solely 
the lowest-order vertices from the effective chiral Lagrangian. The resulting parameter-free expressions for the long- 
range and intermediate-range contributions of types (a), and (b), (c), respectively, can be found in Ref. [28jj, see also 
[29l | where the leading one-loop corrections to the longest-range two-pion exchange terms of type (a) are calculated 
within the infrared-regularized version of chiral perturbation theory. N 3 LO contributions to shorter-range diagrams 
of types (d) and (e) as well as the leading relativistic corrections are given in Ref. [3(|, see also [3l[ for a related 
work. Notice that these shorter-range terms are driven by the leading four-nucleon contact operators which also 
contribute to nucleon-nucleon S-wave scattering. Thus, they do not involve any unknown LECs. Finally, there are 
no corrections to the purely short-range topology of type (f) at this order. An interesting feature of the N 3 LO 3NF 
corrections is their rather rich isospin-spin-momcntum structure emerging primarily from the ring topology (c) in 
Fig. 1. This is in contrast with the quite restricted operator structure of the N 2 LO 3NF. The impact of these novel 
3NF terms on nucleon-deuteron scattering and nuclear structure observables is unknown which makes the complete 
N 3 LO calculations even more urgent, especially in connection with the already mentioned unsolved puzzles. Numerical 
implementation of the new terms in the 3NF at N 3 LO requires their partial wave decomposition which is a nontrivial 
task. In Ref. f32j , a novel method to perform partial- wave decomposition of any type of the 3NF by carrying out 
five-dimensional angular integrations numerically was introduced. This approach is quite general in the sense that it 
can be applied to any type of 3NF but requires substantial computational resources. The partial wave decomposition 
of the N 3 LO 3NF using this new technique is in progress, see Ref. [33| for some first (but still incomplete) results. 

Meanwhile, one may ask whether the derived expressions for the 3NF at subleading order in the chiral expansion are 
already converged or, at least, provide a reasonable approximation to the converged result. This applies especially 
to new operator structures emerging from the genuine loop topologies (b) and (c), whose chiral expansion starts at 
N 3 LO rather than N 2 LO. At this order, the resulting contributions still miss physics associated with intermediate 
A(1232) excitations. In the standard chiral EFT formulation based on pions and nucleons as the only explicit degrees 
of freedom, all effects of the A (and heavier resonances as well as heavy mesons) are hidden in the (renormalizcd) 
values of certain LECs starting from the subleading effective Lagrangian. The major part of the A contributions to 
the nuclear forces is known to be well represented in terms of resonance saturation of the LECs 02,3,4 accompanying 
the subleading ttttNN vertices [2l|, 134143a ] (see, however, the last two references for examples of the A-contributions 
that go beyond the saturation of 02,3,4). The values of these LECs are known to be largely driven by the A and 
appear to be large in magnitude. As a consequence, one observes a rather unnatural convergence pattern in the chiral 

with by far the strongest contribution emerging 

from the formally subleading triangle diagram proportional to C3 



37] 



The (formally) leading contribution to V^ N 
does not provide a good approximation to the potential so that one needs to go to higher orders in the chiral expansion 
and/or include the A- isobar as an explicit degree of freedom. One expects similar conver gen ce pattern for the chiral 
expansion of the 27T-17T exchange and ring 3NF topologies, see also the discussion in Ref. J38|. For the ring topology, 
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this expectation is in line with the phenomenological study of Ref. [39( . All this suggests that one should not truncate 
the chiral expansion of the 3NF at N 3 LO but rather go to (at least) N 4 LO in the standard A-less EFT approach 
and/or include the A-isobar as an explicit degree of freedom. In the latter case, first contributions of the A to the 
27T-17T exchange and ring 3NF topologies would appear already at N 3 LO. It should be understood that the strategies 
outlined above are, to some extent, complementary to each other. This is because N 4 LO 3NF corrections in the A-less 
theory only take into account (some) effects due to single A-excitation but not due to double and triple A-excitations 
which appear first at N 5 LO and N 6 LO, respectively. While these effects are included at N 3 LO in the A-full approach, 
N 4 LO contributions not related to A-excitations are certainly not. We further emphasize that in both cases a number 
of unknown LECs will appear. It remains to be seen which strategy will turn out to be most efficient in practical 
terms. 

In our recent work [4(| we already made a first step in this direction and worked out N 4 LO corrections to the longest- 
range 27r-exchange topology in the delta-less approach. Apart from relativistic corrections (which in our power 
counting scheme appear at N 3 LO but turn out to vanish at N 4 LO), the general form of the 27r-exchange 3NF can be 
parametrized in terms of two scalar functions A(q2) and B{q2) which depend on the momentum transfer q2 = |^| of, 
say, the second nuclcon. In spite of this simple structure, this topology turns out to be most challenging to calculate. 
The pion-nucleon scattering amplitude enters here at the sub-leading one-loop order so that the N 4 LO correction 
depends not only on the pion decay constant F v and the pion-nucleon coupling g^NN but also on 13 independent 

(2} — — — — — — 

(linear combinations of the) LECs from higher-order effective Lagrangians: 01,2,3,4 from d\ +d,2, e?3, d$, du — dig 

from C!~L and §14,15,16,17,18 from The explicit form of the heavy-baryon pion-nucleon effective Lagrangians £^ 

of chiral dimension n needed in the derivation can be found in [4fj| while the complete pion-nucleon Lagrangian 
is constructed in Ref. [4l|. In order to determine these LECs we re-analyzed pion-nucleon scattering at subleading 
one-loop order employing exactly the same power counting scheme as in the derivation of the nuclear forces. We used 
the available partial wave analyses of the pion-nucleon scattering data to determine all relevant LECs. With all LECs 
being fixed from pion-nucleon scattering as discussed above, we found a good (reasonable) convergence of the chiral 
expansion for the functions A(q2) ($(92))- This is to be expected given that effects of the A-isobar are, to a large 
extent, accounted for already in the leading contribution to A(q2) and B(q2) at N 2 LO through resonance saturation of 
the LECs C3.4. As pointed out above, this situation is different for the 2-k-I-k exchange and ring 3NF topologies, whose 
leading contributions at N 3 LO completely miss effects of the A-isobar which lets one expect large N 4 LO corrections. 

In the present work we calculate the intermediate-range contributions to the 3NF at N 4 LO, namely the ones corre- 
sponding to diagrams (b) and (c) in Fig. [TJ and analyze in detail convergence of the chiral expansion for long-range 
tail of the 3NF by comparing the coordinate-space potentials associated with individual isospin-spin-position struc- 
tures. In order to carry out such a comparison in a meaningful way, we worked out the most general structure of 
a local isospin-invariant 3NF both in momentum and configuration spaces and defined the minimal sets of linearly 
independent operators. Our paper is organized as follows. In section [TT] we carry out Fourier transformation of the 
momentum-space expressions for the 27r-exchange 3NF of Ref. Sections Mil and llVl are devoted to the calculation 
of the N 4 LO corrections to the 27r-l7r-exchange and ring topologies, respectively. For 2 ir-lir- exchange contributions 
we provide results both in momentum and coordinate spaces. For the ring topology we give compact expressions in 
coordinate space while the rather lengthy expressions in momentum space are delegated to appendix The most 
general operator structure of a local 3NF is worked out in section [V] where we also define the basis of 22 isospin-spin- 
momcntum operators. We use corresponding coordinate-space version of this basis when discussing numerical results 
for various potentials in section IVII in connection with convergence of the chiral expansion. The findings of our work 
are briefly summarized in section IYIII 



II. TWO-PION-EXCHANGE 3NF IN CONFIGURATION SPACE 



The 27r-exchange topology (a) generates the longest-range contribution to the 3NF. In the isospin and static limits, 
i.e. the limit of infinitely heavy nucleons, its general structure in momentum space has the following form (modulo 
terms of a shorter range corresponding to other topologies): 

VWGfi, qs) = [ g 2 + A j2] +M 2 } ( Tl " r3 - 4 fe) + Tl x Ta ' T2 & x 23 -^B(q 2 )j , (2.1) 

where stays for the pion mass, &i denote the Pauli spin matrices for the nucleon i and ql = Pi' —pi, with pi ' and pi 
being the final and initial momenta of the nucleon i. Here and in what follows, we use the notation: qi = \qi\. Notice 
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that the momentum transfers arc not independent and related to each other via the condition qi + q 2 + q 3 = 0. The 
quantities A(q 2 ) and B(q 2 ) in Eq. (|2.ip are scalar functions of the momentum transfer q 2 of the second nucleon whose 
explicit form is determined by means of the chiral expansion, i.e. the expansion in powers of the soft scale Q ~ M T . 
Unless stated otherwise, the expressions for the 3NF results are always given for a particular choice of the nucleon 
labels. The complete result can then be found by taking into account all possible permutations of the nucleons 

V^ 11 = ^ 3N + 5 permutations . (2.2) 

The explicit expressions for the functions A(q2) and B(q2) at first three nonvanishing orders in the chiral expansion, 
i.e. N 2 LO [Q 3 ], N 3 LO [Q 4 ] and N 4 LO [Q 5 ] 1 are given in Ref. 0. The functions A(q 2 ) and B(q 2 ) resulting at different 
orders in the chiral expansion are plotted versus the values of q 2 in Fig. 5 of that work. While in the case of the 
27r-exchange topology it is possible to address the convergence of the chiral expansion in momentum space thanks 
to the particularly simple parametrization in Eq. (|2.1|) . this is generally not possible for the more complicated cases 
of the 27T-17T exchange and ring diagrams. This is because there is, in general, no easy way to separate the truly 
long-range components, which are unambiguously predicted in terms of the chiral expansion, from scheme-dependent 
short-range contributions. Such a separation is naturally achieved by looking at the corresponding coordinate-space 
potentials at sufficiently large distances. It is, therefore, advantageous and, in fact, also quite natural to switch to 
coordinate space in order to study the convergence of the chiral expansion for nuclear forces. 

We define the coordinate space representation of a static 3NF by means of the Fourier-transform 

V 3N (f 12 , r 32 ) = f^^s J*** ***** W> &)■ (2-3) 
For the two-pion-exchange contribution, we obtain from Eq. (|2.1j) 

VW(rL2, ?"32 ) = -<?l • Vl2 & 3 ' V32 \ Tl ■ T3 -4(r*i 2 ,r32) - Ti X T 3 ■ T 2 V12 X V32 ■ O 2 B(fl2,f32)J , (2.4) 

where = r\ — fj denotes the distance between the nucleons i and j. The differential operators V,j are defined in 
terms of dimensionless variables Xij = ; the functions A and B are given by 



d 3 qi d 3 q 3 , s , 1 1 



(271)3(271)3- ql + MUl + Ml 

d 3 qi d 3 q 3 AS .. ? ,„ AS „. ? „„ 1 1 



A(q 2 ), 



B(q 2 ). (2.5) 



(2tt)3 (2tt)3 " ql + MUl + Ml 

corresponding to A^ 3 \q 2 ) and B^\q 2 ) from Eq. (3.5) of Ref. [4(J are given by 



1287r 2 Fi 



(2c 3 - 4ci - c 3 (Vi2 + V 32 ) 2 ) t/i(a;i2)C/i(x32), 



= ^^V x {x X2 )XJ x {x 32 \ (2.6) 
where §a denotes the nucleon axial vector coupling and the Yukawa function U\ is defined as 



-4(n.2,r 32 ) 
B{r 12 ,f 32 ) 
The N 2 LO expressions for A and B 
A {3) (r 12 ,r 32 ) 
B {3 \r 12 ,r 32 ) 



4tt f d 3 q e^/^ _ e~ x 
~IUJ (2tt)3 g 2 + M 2 ~ — ■ 



(2.7) 



Here and in what follows, the superscripts of A, B, A, B as well as other functions parametrizing 3NF matrix elements 
refer to the chiral dimension, i.e. to the associated power of the soft scale Q. 



1 Notice that the overall "chiral dimension" is a matter of convention. In the context of nuclear chiral EFT, one usually uses the convention 
in which the leading-order (LO) one-pion exchange nucleon- nucleon potential is assigned the chiral dimension Q°. 



The first corrections to A and B emerge from Fourier-transforming the expressions A^\q 2 ) and B^\q 2 ) given in 
Eq. (3.4) of Ref. 0. We obtain 

-A (4) (ri2,r 32 ) = 409^6 { [ M + l) - 2 (<£ + l) (V 12 + V 32 ) 2 ] C/ 1 (x 12 )t/ 1 (x 32 ) 

+ i- (2 - 5(V 12 + V 32 ) 2 + 2(Vi 2 + V 32 ) 4 ) / d 3 x [7i(|2 12 + £|) ETiflSaa + f|)}, 

^ (4) (r 12 ,r 32 ) = - i ||-^{(2 5 i + !)[/, (x 12 )t/ l( x 32 ) 

+ ^ ( 4 - (^12 + V 32 ) 2 ) / d 3 x [7i(|2 12 + S\) W^x) l/i flfia + 2|)}. (2.8) 

The profile function Wi is given in terms of the Fourier-transform of the loop function A(q) appearing in A^(q 2 ) 
and B^(q 2 ): 

4tt f d 3 q m(Ui! e~ 2x Af _^ 1 q 



WAx) = —5 / — e lq - x ' M - A(q) = =- , with A(q) = — arctan — 1 — . (2.9) 

U ' Ml J (2tt) 3 w 2x 2 ' w 2q 2M W V ; 

To give the coordinate space expressions for N 4 LO contributions we need to Fourier-transform another loop function, 
namely 



l(9)= VF±M log v?+M±«, (2 ,o) 

which enters the expressions for A (5) {q 2 ) and S (5) (g 2 ) in Eq. (3.14) of Ref. 0. This can be most easily achieved by 
using the spectral representation of L given by 

f°° n 2 1 

L (1) = 1 + / ^^2^VM 2 -4M2. (2.11) 
The Fourier-transform of the square-integrable part of L is given by 

a r ,73 „ /*oo 1 -1 -1 />oo — x fj. 

With these preparations, we obtain the following result for the Fourier-transform of A^ 5 \q 2 ) and B^ 5 \q 2 ): 
i (5) (ri 2 ,f 32 ) = 73 ^f F6 [ - (Vi 2 + V 32 ) 2 (F 2 (23047r 2 5A (4e 14 + 2e 19 - e 22 - e 36 ) - 2304rr 2 J 18 c 3 ) 

+ 5 a(144ci - 53c 2 - 90c 3 )) + F 2 (4608tt 2 J 18 (2 Cl - c 3 ) + 4608 7 r 2 < M (2e 1 4 + 2e 19 - e 36 - 4e 38 )) 
+ g A (72 (647r 2 r 3 + l) d - 24c 2 - 36c 3 ) 

+ (Vi 2 + V 32 ) 4 (23047r 2 ei 4 F 2 ff J 4 - 2<m(5c 2 + 18c 3 )) ] f/i(o;i 2 )f/i(a; 32 ) 
9\Ml 



12288tt 4 F6 



1 - 2(Via + V 32 ) 2 4(6ci - c 2 - 3c 3 ) - (V12 + V 32 ) 2 (-c 2 - 6c 3 ) [^12)^1(3:32) 



49152 



-f^ (l - 2(Vi2 + V 32 ) 2 ) (4(6 Cl - c 2 - 3c 3 ) 



- (Via + V 32 ) 2 (-c 2 - 6c 3 )) (Via + V 32 ) 2 / d 3 x [/ifljfo + x\) Vx{x) J7i(|x 32 + 

6 (5) (ri 2 , r 32 ) = - 36864^4^ 6 [ ^ ( 1152 ^ d ~ 18C4 ~ H527rV(2ei 7 + 2e 2 i - e 37 )) + 108^c 4 + 24g A c 4 

- (V12 + V 32 ) 2 (5 5A c 4 - U52n 2 e 17 F 2 g A ) ^(x^U^) 

+ ||^( 4 -^ 2 +^ 2 ) 2 ) c/i ^ 2 ) c/i ^ 2 ) 



24576tt 5 F6 



(4 - (Via + V 32 ) 2 ) (Via + V 32 ) 2 / d 3 x E/i(|#i2 + £|) Vi(a:) t/i(|f 32 + f|). (2.13) 
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It remains to emphasize that while the momentum space representation of the functions A and B depends on just 
one variable q 2 , the coordinate-space functions A and B depend on three scalar arguments. We will discuss the 
convergence of the chiral expansion for the coordinate space potentials in section IVII 



III. TWO-PION-ONE-PION EXCHANGE 3NF AT N 4 LO 



We now turn to the 2tt-\t: exchange topology. In contrast to the longest-range 2-k exchange topology discussed in 
the previous section, its chiral expansion starts at N 3 LO. At this order, one has to evaluate all one-loop diagrams 
made out of the lowest-order pion-nucleon vertices. This was achieved in Ref. (28|, see Eqs. (2.16)-(2.23) of that work. 
As pointed out in Ref. [io| . the decomposition of momentum-space 3NF expressions according to the type of the 
topology is not unique as e.g. some parts of the 27r exchange contributions can be reshuffled into 2-k-I-k exchange and 
shorter-range terms by canceling pion propagators with the corresponding expressions in the numerator. In Ref. (40| 
we introduced a "minimal" parametrization of the 2tt exchange 3NF which corresponds to Eq. (|2.ip and which is 
adopted here and in what follows. 

The structure of two-pion-one-pion exchange contributions up to N 4 LO in the chiral expansion has the form 



^3 • g3 

- 2 1 Ml 



Ti ■ t 3 [a 2 ■ qi qi • <?3 F 1 (q 1 ) + a 2 ■ q\ F 2 {qi) + a 2 ■ q 3 F 3 (qi)] + t 2 -t 3 [Si ■ q\ q\ ■ q 3 F^qi) 

q 3 + 1V1- L 

+ c?i ■ q 3 F 5 (qi) +a 2 ■ qi q\ ■ q 3 F 6 (qi) + 5 2 ■ qi i*V(gi) + a 2 ■ q 3 q\ ■ q 3 Fg(qi) + a 2 ■ q 3 F 9 (qi)] 

+ Ti x t 2 ■ 1-3 [<?i x a 2 ■ qi {qi ■ q 3 F w (qi) + Fu(qi)) + q\ X q 3 ■ o x q\ ■ a 2 F 12 (qi)] , (3.14) 

where i 5 !... 12(91) are scalar functions to be calculated. Notice that we use a slightly different notation compared to 
our early paper [28l| . which is now also valid at N 4 LO. First non- vanishing contributions to the structure functions Ei 
are generated at N 3 LO by diagrams shown in Fig. 3 of that work. Adjusting the expressions obtained in Ref. [28| to 
our new notation and taking into account terms induced by reshuffling the 27r exchange contributions as explained 
above we obtain the following results for the functions Fj(gi): 
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^ 4)(9l) = lSF% A{qiH2M " + ql) ' 

*s W («0 = -^^A{ qi ){{S 9 \~A)Ml + {i 9 \-l)ql), 

F A ^ q^~-^Fl A{q ^ 

F^\ qi ) = Ft\ qi ) = F^\q 1 )=F$( qi ) = F£\q 1 ) = Q, 



F ^ {qi) = _ 5I2^ A((?l) (4M ' + ql) ■ (3 - 15) 



Notice that we give here only non-polynomial parts as the polynomial ones simply lead to shifts of the low-energy 
constants D and E from N 2 LO three-body force. 

First corrections to these results emerge at N 4 LO from diagrams shown in Fig. [21 which involve a single insertion of 
Ci-vertices from the sublea sing pio n-nucleon Lagrangian. Evaluating the irreducible contributions of these diagrams 
following the lines of Refs. [28l . |3G ] and keeping only terms non-polynomial in qi we obtain the following expressions: 

F ^ iqi) = - ^F^l{AM i + ql) L{qi) ( 8 - $ M * + 2 ^ + $ Mlql ~ & - ^ ^ 
{\-lg\)g\c A Ml 
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(4) 



(5) 



(6) 



(7) 



(8) 



(9) 



O V- 



(10) 



FIG. 2: Two-pion-one-pion exchange 3N diagrams at N 4 LO. Solid dots and filled circles denote the leading and subleading 
pion-nucleon vertices, respectively. Graphs resulting from the interchange of the nucleon lines and/or applying time reversal 
operation are not shown. For remaining notation see Fig. Q] 



4V) 



*i ( o%i) 



48tt 2 F6 (4A/2 + q\ 



-Mqi) (4 (4g 2 A - l) + (l7g 2 A - 5) M 2 q\ + (4g 2 A - l) qf] 



16tt 2 F6 



L{qi), 



g A M^{6 Cl + c 2 - 3c 3 ) 

967T 2 F6g 2 

192 ^4^(4M2 + g 2 ) H 8 ^ 4 + C2 (-8Jl# + 2M 2 q 2 + qf) + 12c 3 A/ 2 (2M 2 + q 2 )) 
- 1Q 2^FS L{qi) ( 24ciM - " C2 ( 4M * + «?) - 6c3 ( 2M - + . 



128tt 2 F6 (4A/ 2 + g 2 ) 
c 3 (80M* + 68M 2 g 2 + 13gJ 



- 32 Cl M 2 (3M 2 + qf) + c 2 (16M* + 16M 2 q 2 + 3<tf) 



4 2 %i) = 



9A c ^ L (qi) 

64tt 2 F6 



(3.16) 



It is straightforward to transform these results into coordinate space. The general structure corresponding to the 
momentum-space expression in Eq. (|3.14[) has the form: 

V27r-l7r(n.2,r32) = 0*3 ' V 32 Tl • T 3 ^0*2 ' Via Vl2 ■ V 32 A (x 12 ) - CT 2 ■ Vl 2 F 2 (x 12 ) - <7 2 ■ V 32 F 3 (x\ 



+ T 2 ■ T 3 ^CTl • Vl2 V12 • V 32 F 4 (xi 2 ) - G 1 • V 32 F 5 (.Ti 2 ) + CT 2 • V12 V i2 • V 32 F 6 (xi 2 ) 
- a 2 ■ V12 -F 7 (a;i 2 ) + a 2 ■ V 32 V 12 • V 32 F s (xi 2 ) - ct 2 • V 32 F g (x 12 )j 
+ ti x r 2 ■ r 3 ^<ti x <r 2 • V12 (Via ■ V 32 Fio(^i2) - ^11(^12)) 

+ V12 x V 32 ■a 1 V12 -0*2^12(2:12))] ^1(^32) • (3.17) 
In order to calculate the functions i*j it is convenient to employ the spectral representations of the function L(q), see 
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Eq. (|2.1ip . and A{q) given by 



1 Z" 00 1 

A{q) = 2 / dM ^T 

1 J2M W M + 



Following Ref. [28|], we define the profile function Ws(x) via 



47T 



d 3 g 



7 (2tt) 3 



where 



Jq-x/M^ 



Ei{x) 



M AM 2 

- ^(«) 



e~ l dt 



2Ei(-2x) 



t 



(3.18) 



(3.19) 



(3.20) 



Fourier transform of terms in Eq. (|3.16[) involving the function L{q) can be expressed using the profile functions V\(x) 
from Eq. (|2.12p and V2(x) which is defined according to 



V 2 (x) = 



dfj,- 



-X fl 



M 2 \/m 2 ~ 4 

With these definitions, the N 3 LO contributions to the Fi-functions are given by 



(3.21) 



5(4 



F. 



p(4 



F 



(4 



(£12) 

(»12) 

O12) 

(»12) 

Ol2) 
O12) 

(^12) 

O12) 



gW 

4096tt 3 F6 

2048tt 3 F6 

4096tt 3 F6 1 

2048tt 3 F6 

2048 ir 3 F$ 
V12) 

2048tt 3 F6 



(2g 2 A U 1 (2x 12 ) - (g 2 A + l) Wi(xia) + (2^ - l) W 3 (x 12 )) 
(V 2 2 -2)^i(x 12 ), 
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(3V 2 2 -8) si + ^VF!^), 



W!{X 12 ), 



V 2 2 VF!(a; 12 ), 



ffW = F 9 W(x 12 ) = F$(x 12 ) = f£\x 12 ) = 0, 



2) Wi{x 12 ), 



Ami 

8192tt 3 F6 



(V 2 2 -4) ^1(112), 



(3.22) 



while for the N 4 LO contributions we obtain the following results: 



Jgf, [(V? 2 - 4) (-V 2 2 + (V 2 2 + 10) £ + 2) E/^) 



12 + iov^ 2 



F^(x 12 ) 



2 (-V} a + 2V 2 2 + (V 
15367^*5 1 12 12 ( 



32) g 2 A + 8)V 2 (x 12 )] 



12 - 17V^ 2 



16) 5i"4) 



F 4 (5) (x 12 ) = 



((V? 2 - 4) t/i(2x 12 ) - 2V 2 2 F 2 (x 12 )) 
fV 2 2 Vi(x 12 ), 



gAC 4 M| 2 



2567r 4 F6 
glc 4 M| ^ 



2567T 4 F6 
,4 



V? 2 yi(x 12 ), 



61 g 4 ^ f6 [2 (48ci + (-Vt a + 2V 2 2 + 8) c 2 + 12 (V 2 2 - 2) c 3 ) y 2 (x 12 ) 
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(V? 2 - 4) ((V? a - 2) cs, - 12c 3 ) C/i(2 a ; 12 )] , 



F 7 (5) (x 12 ) = 



AMI ^ 



V? 2 (24ci + (V? 2 - 4) c 2 + 6 (V? 2 - 2) c 3 ) 7i(i 12 ), 



3072tt 4 F6 

F 9 5 \xii) = JnaJj* (32 (V? 2 - 3) Cl + (3Vf 2 - 16V? 2 + 16) c 2 + (l3Vf 2 - 68V? 2 



4096^^ 

x ((V? 2 -4)t/ 1 (2x 12 )-2V? 2 y 2 (x 12 )) 



C3 



g A c 4 A^ 2 



(3.23) 



IV. RING DIAGRAMS AT N 4 LO 



Finally, we consider the ring topology. The leading contributions emerge at N 3 LO from diagrams shown in Fig. 4 of 
Ref. [281 ] . As explained in that paper, only diagrams proportional to g\ and g\ generate nonvanishing 3NFs: 



y{4) = y(4),9i y{4),A 
ring ring ' ring 



(4.24) 



Evaluating the corresponding loop integrals in momentum space we obtained complicated expressions involving three- 
point function which are given explicitly in appendix of Ref. [28j. 2 The results in coordinate space are much more 
compact and have the form: 



Kin, A (ri2,r 32 ) 



QA \ 6 f d 3 L cPIo d?l 3 J? 2 if f) ;? a V 

2F„ J J (2tt)3 (2tt)3 (2tt)3 [If + M%] \l\ + M%] 2 [|§ + M%] 

Ami 



4096 71-3^6 



4ti ■ t 2 V 23 x V12 ■ ct 2 V 23 x V 3 i ■ o 3 V 3 i ■ Vl2 



2ti • r 3 V 23 • V 3 i V 23 • V12 V 3 i • V12 + ri x r 2 • t 3 V 23 x V12 ■ <? 2 V 23 • V 3 i V 3 i • V i5 



- 3V 3 i x V12 • <?i V 23 x V 3 i • 3 3 V 23 • V12 



C/i(a;2 3 ) U 2 (x 31 ) Ui{xi 2 ) , 
2n ■ t 2 (V 23 • V 3 i V 3i ■ V12 - V 3i x V12 • cri V 23 x V 31 . g 3 ) 

Ui{x 23 ) Ui(x 31 ) f7i(xi 2 ), 



2048 7T 3 F6 

+ Ti x r 2 • r 3 V 3i x V12 • CTl V 23 • V 3 1 



(4.26) 



where the derivatives should be evaluated as if the variables x\ 2 , x 23 and x 3 \ were independent 3 and the numerator 
v in the first line is given by 

v = -8ti ■ t 2 h x l 3 ■ ct 2 h x l 2 ■ S 3 h-h - 4ti • r 3 ii • l 2 h ■ h h ■ h + 2ti x t 2 ■ t 3 h x l 3 ■ o 2 h ■ l 2 l 2 ■ l 3 

+ 6f 2 x l 3 ■ g i h x h ■ ct 3 h ■ l 3 . (4.27) 



At N 4 LO, one only needs to evaluate the contributions of the four diagrams shown in Fig. |3l 



V, 



(5) = y{S),g\ y(5),g 2 A y(,5),g A 
ring ' ring ' ring 



(4.28) 



2 We emphasize that several symmetry factors are missing in Eq. (Al) of that work. The corrected equation has the form: 

1 1 t _, 

-a\ ■ 530-3 ■ gi Kg 



0*1 ■ CT2 T2 ■ T3 R% + CTl ■ lji(T2 ■ qi T% • T3 R2 + S\ ■ q\U2 ■ 93 r% ■ T3 R3 + <?i ■ 93^2 • 91 T"2 • T3 J?4 

1 

— f 

2 



<?1 ■ <?3<52 • 93 T2 • X3 R 5 + -Tl ■ T3 /? 6 + CTi ■ (J1CT3 ■ (fi fly + -CTi ■ q\(f a ■ q S 



-cti ■ e?3 J?io + -gi ■ <?*3 x ct 2 ti ■ t 2 x t 3 fin. 



(4.25) 



Clearly, the relative distances r*i2, r*23 and r3i are related via r*i2 + r23 + r3i = 0. 
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(1) (2) (3) (4) 

FIG. 3: Ring diagrams at N 4 LO. For notation see Figs.CQE] 

We were again able to obtain fairly compact expressions in coordinate space, which, however, involve now a single 
scalar integral over the mass of the exchanged particles: 



•ring 



Ami 



ds 



2V 12 • V 



12 ' V 23 1 V 23 



V 3 i (12cict 2 • a 3 - 4c 2 s 2 (T2 • (T3T2 • t 3 + 6c 2 s 2 



10247r 4 i^ 

+ V12 • V 3 i((Ti • <? 3 (-3c 3 + C4X1 • t 2 + c 4 r 2 • r 3 ) + 2(<t 2 • <? 3 (-3c 3 + c 4 Ti • t 2 + c 4 Ti • t 3 ) 

+ c 3 (2tx ■t 2 + t 1 - r 3 ))) + 2V 12 • 0W31 • <7 2 (3c 3 - c 4 t 1 ■ t 3 - c 4 r 2 • r 3 ) - 4c 3 s 2 <7 2 • a z r 2 ■ t 3 

+ 6c 3 s 2 + c 4 V 12 • V 31 x d X Ti ■ t 2 x t 3 ^ - 2 ^2Vi2 • <?iV 31 • a 2 (3c 4 - s 2 (c 2 + c 3 )ti • r 2 ) 

+ fv 3 i ■ CT1V31 ■ (?3 - (s 2 + l) cti ■ (T 3 ) (s 2 (c 2 + c 3 )xi ■ t 3 - 3c 4 ) + V12 • V 3 i(4ciTi ■ r 2 

+ V 23 ■ CT1V31 ■ a 3 (-3c 3 + C4T1 • -r 2 + c 4 r 2 ■ t 3 ))) - (s 2 + l) V 12 ■ <t 3 V 23 ■ <?i(3c 3 - c 4 r 4 • r 2 

- c 4 r 2 • t 3 )) + V12 • V31 (8 (v 12 • CT2V23 • B\ (s 2 (c 2 + c 3 )ri • t 2 - 3ci) 
+ V 23 ■ CT1V31 ■ a 3 (s 2 (c 2 + c 3 )ti ■ r 3 - 3ci) + (s 2 + l) a 2 ■ a 3 (s 2 (c 2 + c 3 )t 2 • r 3 - 3ci) 

- V23 • V31 \Aai • <r 3 (s 2 (c 2 + c 3 )ri • r 3 - 3c 4 ) + 8ciTi ■ r 3 - 6c 2 s 2 + 4Vi2 ■ CT2V23 ■ e?i 
x (-3c 3 + c 4 x 4 • r 3 + c 4 r 2 • t 3 ) - 6c 3 s 2 + c 4 Vi 2 ■ V23 x <t 2 ti • t 2 x t 3 



4 ( (s 2 + 1) V 23 • v\ (V31 ■ o 2 (6ci - 2s 2 (c 2 + c 3 )ti • r 2 + V 23 ■ V 3 i(-3c 3 + c 4 ti ■ t 3 + c 4 r 2 • r 3 ) 



V12 ■ o 3 (3ci - s 2 (c 2 + c 3 )ti ■ t 3 ) J + V12 ■ (T1V12 ■ (T2V23 ■ V31 (6c 4 - 2s 2 (c 2 + c 3 )ri ■ r 2 
V23 • V 3 i(-3c 3 + c 4 ti • t 3 + c 4 r 2 ■ T 3 )j J - 2(Vi2 • V 23 ) 2 (V31 • CT1V31 • <? 3 - (s 2 + l) &x ■ a 3 



(3C 3 - C 4 T 4 ■ T 2 - C 4 T 2 ■ T 3 ) + 4 (,S 2 + 1) (Vl2 ■ V 3 l) 2 (T2 ■ <7 3 (3c 3 - C 4 Tl ■ T 2 - C 4 Tl • T 3 ) 



V. 



x C/r(a;i2)C/r(x23)C/2 S ^3i) 



(4.29) 



(5),si 



1024^^6 ^ 



ds 



8C1V12 • V23T2 ■ T3 + 8ClV 23 ■ V31T2 ■ t 3 + 4c 2 s V12 • V23T2 • T 3 



V 



(5),9° 4 



ring 



+ 4c 2 s 2 V 2 3 • V31T2 • t 3 + V12 • V31 (-4c 3 V 12 • V23X2 • t 3 - 4c 3 V 23 • V 3 ir 2 • r 3 

+ C 4 Vl2 ■ V 2 3 X <7 2 Tl ■ T 2 X T 3 + 4C 4 V 2 3 ' OlV31 • 3 (j X -T 2 +T 2 - T 3 ) 

+ 2c 4 (s 2 + l) <7 2 • 53(1-1 -t 2 + ti ■ t 3 )) + 4c 3 s 2 V 12 ■ V23T2 • t 3 + 4c 3 s 2 V 2 3 • V 3 ir 2 • t 3 

- 2c 4 Vi2 ■ V 2 3V 3 i ■ CT1V31 ■ CT3T1 ■ t 2 - 2c 4 Vi2 ■ V23V31 • (T1V31 • a 3 r 2 ■ r 3 

- 2c 4 Vi2 • V23V23 ■ V31CT2 ' CT3T1 ■ T 2 - 2c 4 Vi2 ■ V23V23 • V 3 i<72 • S 3 Ti ■ T 3 

+ C 4 V 2 3 • V 3 lVl2 • V 23 X 2 T\ ■ T 2 X T 3 - 2C 4 (s 2 + l) V12 • CT 3 V 2 3 • <?l(Tl • T 2 + T 2 ■ T 3 ) 

x Ut(x 12 )Ut(x 23 )Ut(x 31 ), 



ds 



(4.30) 



10247r 4 i^ 



7T J —OO 



4r 2 • r 3 ( 2ci + s (c 2 + c 3 ) - C3V12 ■ V31 + C4V12 • V23 x a 2 Ti • r 2 x r 3 
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x U?(x 12 )UZ(x 23 M(x 31 ), (4.31) 

where 



The expressions in momentum space are rather lengthy and can be found in appendix El 



V. GENERAL OPERATOR STRUCTURE OF A LOCAL THREE-NUCLEON FORCE 



As already emphazized in the introduction, having derived explicit expressions for the long-range part of the 3NF at 
the three first orders in the chiral expansion, it is interesting to test convergence in coordinate space. One generally 
expects for the chiral expansion of nuclear potentials to converge at distances of the order of or larger than r ~ M~ x . 
In order to analyze the convergence of the chiral expansion for three-nucleon potentials in a meaningful way, we 
first need to define a basis in the space of isospin-spin-position or, equivalcntly, isospin-spin-momcntum three-nucleon 
operators. Thus, we need to work out the most general structure of the three-nucleon force. To the best of our 
knowledge, this task has not been accomplished yet, see however Ref. [42j, where the most general isospin structure 
of the 3NF is given. 

Given that a general 3NF depends, in the center of mass system, on four independent momenta in addition to the 
spin and isospin Pauli matrices, its structure is obviously rather rich. Fortunately, even at such a high order in the 
chiral expansion as N 4 LO, the most complicated part of the three-nucleon force (before antisymmctrization) is still 
local. For the long-range part, the only non-localities in the power counting scheme we adopt arise from the leading 
rclativistic corrections to the 27r exchange diagrams discussed in Ref. [30(. We, therefore, restrict ourselves here to 
the most general structure of a local 3NF. We, furthermore, require in the following that the 3NF is invariant under 
parity, time-reversal and isospin transformations. 

Every operator appearing in the 3NF can be written as a linear combination of spin-momentum terms multiplied 
with isospin structures. We remind the reader that according to the standard convention, the expressions for nuclear 
forces are to be understood as matrix elements with respect to momenta and operators in the spin and isospin spaces. 
The building blocks for the spin-momentum structures are 

<?i, ct 2 , <?3, q\, q 3 , (5.33) 

where cr, are the Pauli spin matrices of the nucleons while q\, q 3 denote the two independent relative momenta. 4 
These building blocks have to be contracted with the tensors S a b and e a b c to build scalar operators. We have the 
following symmetry constraints: 

• Parity invariance of the force allows only for spin-momentum structures which arc invariant under 

qi -> -qi and q 3 -t -q 3 . 

• Time-reversal invariance implies that only those structures contribute which are invariant under 

fii -> —Bi , qi ->• Qi and rf -> — rf , i = 1,2, 3, 



see Eq. (2.47) of Ref. (43 1 



4 The momentum transfer g*2 can be expressed in terms of gi 3 via g*2 = —Qi — Q3- 

5 The invariance under rf — > —rf follows directly from the invariance of the matrix element under (t'\Ti\t) (t\Ti\f) 
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• Isospin conservation requires any structure to be a product of a spin-momentum operator with one of the 
following isospin-structures: 

1, n-T2, T\ ■ T3, r 2 ■ T3 and n • (r 2 x r 3 ). 
In addition to the symmetry constraints, we will also employ the Schouten identity 

to eliminate redundant structures. A general local three-nucleon force can be written in a form 

^ Oi(ai, &2, V3,Ti,t 2 , r 3 , qi, qk) Fjjqi, g3, gi ■ 93) , 

where Oi arc spin-momentum-isospin operators and the scalar structure functions Fi depend only on absolute values 
|gi|j 1 93 1 and on the scalar product q\ ■ q^. The three-nucleon force V^ 11 in Eq. ()2.2j) is obviously invariant under any 
permutation P € S3, with S3 denoting the permutation group: 

Y,POiPF i = Y J O i F i , (5.34) 

i i 

where 

POi(a 1 ,a 2 ,ff3,T 1 ,T 2 ,T 3l q 1 ,q 3 ) = Oj(a P [i] , a P [ 2 ] , op [3] , Tppj , rppj, r P [ 3 ] ,qp[i], <?p[ 3 ]), 

PFi(qi,Q3,qi ■ 93) = ^i(gp[i],gp[3]>g>[i] • g*P[3])- (5-35) 

To understand the behavior of the structure functions under permutations of momenta it is necessary to analyze the 
behavior of the operators Oi under permutations. Since the operator set we consider here is complete, the permuted 
operator POi is just a linear combination of Oj's: 

PO i = J20 j D ji (P), 

3 

where D(P) are some invertible matrices. It is easy to see that the set of matrices D builds a representation of S3. 
Indeed 

P'POi =P'J2 O k D kl {P) = °i D 3k(P') D kt (P) = £ (), !),,■ I" l':. 

k j,k j 

from which immediately follows 

D(P'P) = D(P')D(P). 

Transformations of the structure functions Fi with respect to permutations P of the nucleon labels can now be read 
off from 

J2 OiFi = £ POiPFi = OjD^PFt =E°4Z Dij(P)PFj , 

i i i,j i \ j ) 

from which we obtain the identity 

F i = J2D ij (P)PF j . 

3 

It is advantageous to choose the basis in the space of operators Oi such that the representation matrices D arc 
block-diagonal corresponding to irreducible representations of the group S3. There are three inequivalent irreducible 
representations of S3: 

• The trivial (identity) and antisymmetric (— 1)™(- P ) representations with w(P) = ±1 for even/odd permutations 
are one dimensional. 
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• The third irreducible representation is two dimensional. The representation matrices can e.g. be chosen as 

(5.36) 



ill), Z>((12)) = § ( ^- 2?((13)) =( o i 



Mm = -h ( ^ ) , ©((123)) = "3 ( ) . C (( 132 )) = -5 ( ^ 

where we used the cyclic notation for permutations: 

0[2] = 2, ()[3] = 3, 

(12) [2] = 1, (12)[3] = 3, 

(13) [2] = 2, (13)[3] = 1, 
(23)[2] = 3, (23)[3] = 2, 
(123)[2] = 3, (123)[3] = 1, 



ON 


= 1, 


(12)[1] 


= 2, 


(13)[1] 


= 3, 


(23) [1] 


= 1, 


(123)[1] 


= 2, 


(132)[1] 


= 3, 



(5.37) 



To construct the operators Oi for which D(P) is block-diagonal, we introduce the symmetrizing and antisymmetrizing 
functions: 

= I E PO > A ^°) = \H (-i) w(p) po. 

Pes 3 Pes 3 

It is obvious that 

PS(0) = S{0) and PA{0) = (-l) w(P) A(0) 

for all P E S 3 such that S(0) and A(0) transform under one-dimensional representations. To construct operators 
which transform under two-dimensional irreducible representation we introduce the functions 



G v (°) = \ m (P)PO, i, 3 = 1, 2. 



3 

Pes 3 

It is easy to verify that the resulting operators Gij (0) indeed transform under two-dimensional irreducible represen- 
tation: 

PGaip) = \ J2 Vij(Q)PQO = \J2 m( p ^Q)QO = J2 G kj(o)v hl (p). 

Qes 3 Qes 3 k 



With all the symmetry constraints introduced above, we found that the most general structure of a local 3NF can be 
written in terms of 89 operators 0\, . . . , 089, which transform with respect to permutations according to irreducible 
representations of S3. These 89 operators can be generated from a set of 22 independent operators Gi, - ■ ■ , G22 using 
the functions S, A and Gy defined above. The explicit form of generating operators Gi, ■ ■ ■ , G22 we adopt in this work 
and their relation to the generated operators 0\ , . . . , 089 are given in Table Q] The complete expression for a local 
three-nucleon force in our notation can be written in the symmetric form 

89 

^3N U = ^20 l (a 1 ,a 2 ,a3,T 1 ,T2,T 3 ,q 1 ,q 3 )F l (q 1 ,q 3l q 1 ■ q 3 ). (5.38) 

i=l 

In this representation, the structure functions Fi have simple transformation properties with respect to permutations 6 . 
An alternative way to express the three-nucleon force is given by 

22 

^3N U = '^2Gi(ffi,^2,^3,Ti,T 2 ,T 3 ,qi,q 3 )F l {q 1 ,q 3 ,qi ' ^ + 5 permutations. (5.39) 

i=l 



The only exception is F17 which mixes different contributions from other structure functions. This is due to the use of the Schouten 
identities. 
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Generators Q of 89 independent operators 


S 


A 


G12 


G22 


Gn 


G21 


01 = 1 


Oi 

















02 = Tl ■ T3 


O2 





o 3 


o 4 








03 = ffl • <?3 


O s 





o 6 


o~ 








04 = Tl • T3<7l ■ (73 


8 





o 9 


O10 








05 = T2 ■ T3CT1 ■ <?2 


On 


012 


O13 


O14 


O15 


Ole 


06 = Tl • (t 2 X T 3 )o*l ■ (<T 2 X <7 3 ) 


O17 

















t/7 = Tl ■ (t 2 x t 3 )ct 2 ■ (qi x g 3 ) 


Ois 





O19 


O20 








08 = qi ■ aiqi ■ S3 


On 


022 


023 


o 24 


Ois 


o 26 


09 = gi ■ ?3?3 ■ G\ 


On 





o 28 


o 29 








0io = gi ' <?1<Z3 • <?3 


O30 





o 3 l 


o 3 2 








011 = T2 ■ T 3 gl ■ Siq\ ■ <?2 


O33 


034 


035 


o 36 


037 


o 38 


012 = T2 • T3<fi ■ 9\<fe ■ Si 


O39 


O40 


o 4 l 


o 42 


043 


044 


013 = T2 ■ T3<73 ■ S\q\ ■ <?2 


O45 


o 46 


047 


o 48 


049 


O50 


014 = T2 • T3g*3 ■ <7i<f3 ■ Si 


Obi 


o 5 2 


053 


o 54 


o 55 


o 56 


015 = Tl ■ T3<f2 ' S\q\ ■ Sz 


O57 





O58 


O59 








016 = T2 • T3IJ3 ■ (729*3 ' S3 


O 6 o 


Oei 


6 2 


6 3 


6 4 


o 65 


017 = Tl ■ T 3 gi ■ CTigjj • S3 


o 66 





o 67 


o 68 








018 = Tl ■ (t 2 X T3)Sl ■ S3S1 ■ (qi x <f 3 ) 


O m 





O70 


O71 








019 = Ti ■ (t 2 x t 3 )i73 ■ q\qi ■ (Si X Si) 


O72 


O73 


074 


O75 


o 76 


077 


020 = Tl • (T2 X T 3 )«7i • qiS'2 ■ gi<7 3 ■ (gi x g 3 ) 


o 78 


O79 


O 80 


Osi 


82 


o 83 


021 = Tl ■ [Tl X T 3 )<7i ■ g 2 <73 ■ g*2<72 ■ (<fl X <f 3 ) 


8 4 





8 5 


o 86 








022 = Ti • (t 2 X T 3 )«7i • qi<7 3 ■ g 3 <72 • (gi X <f 3 ) 


87 





8 8 


o 89 









TABLE I: The set of 22 generating operators Qi and their relation to 89 independent operators Oi , . . . , Os9 which parametrize 
the most general structure of a local 3NF. The operators Oi are generated by application of one of the 6 functions 
5,74,011,012,021,022 on the corresponding operator Qj. The 22 operators are constructed to be either totally symmet- 
ric, symmetric under 1 f-> 3 or unsymmetric. 



It is easy to see that the two representations (|5.39|) and (|5.38l) are equivalent. While Eq. (|5.39p can obviously be 
brought into the form of Eq. (|5.38p we now show that the converse is also true. Eq. (|5.38[) can be rewritten in the 
form 



22 



viif 1 = J2{ S (^)M l + A(g i )N l + ]T G jk (gi)L 



J,fe=l 



22 f 2 1 

= E 5>(&) \ M * + J(-ir (p) ^ + E \v jk (p)L% 



(5.40) 



where Mi, iVj, are some of the structure functions Fi(l — 1, . . . , 89) from Eq. (|5.38|) . From the symmetry property 
(|5T34)l we get 



v ™ = I E p '( y ™) 

P'es 3 

22 f 2 \ 

= 1 E f^P'P(Gi)\lP'(M i ) + h-ir^P'(N i )+ £ \v jk {P)P>{ti> jk )\ 

P',P>eS 3 i=l [ j,k=l J 

22 f 2 | 

= E P "X> E { 4 p ~ 1(Mi) + ^(-v^p-Hm + E k D ^{p) p ~\ L %) \ > ( 5 - 41 ) 
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where we made a change of variable P" = P'P in the last line. This equation has the form of Eq. f|5 .39[) with 

K ■= E I ^P- 1 (M i ) + ±(-ir^p-\N i )+ £ ^V jk (P)p-\L) k ) 1 . (5.42) 
Pes 3 [ j,k=i J 



VI. CHIRAL EXPANSION OF THE LONG-RANGE TAIL OF THE 3NF 



With these preparations we are now in the position to address the convergence of the chiral expansion for the long- 
range tail of the 3NF. It is clear that all arguments of the previous section can also be applied to operators in 
coordinate space. Here and in what follows, we use the following basis of 22 operators: 



0i 


= 


1. 
















02 


= 


Tl ■ 


1"3 , 














g 3 


= 


#1 ■ 


fa , 














04 


= 


Tl ■ 


T 3 fa 


' ^3 , 












Gs 


= 


T 2 ' 


1~3 fa 


■ ^2 , 












06 


= 


Tl ■ 


(r 2 X 


T 3 ) Oy ■ 


(o 2 x 


fa) , 










= 


Tl ■ 


(T2 X 


T 3 ) O 2 ■ 


(ri 2 x 


f2 3 ) , 








Gs 




f 23 


■ (71 f 23 


■ o 3 , 












09 




f 23 


■ fa f 12 


■ Ol , 












Gw 




?23 


• #1 fl2 


■ o 3 , 












Gn 




T2 ■ 


T"3 ^23 


■ Oi f 2 3 


• o 2 , 










Gl2 




T2 ' 


T 3 f 2 3 


■ Ol fu 


■ o 2 , 










013 




T2 ' 


T 3 fi2 


■ Ol r 23 


■ o 2 , 










014 




T2 ' 


T"3 fi2 


■ Oi r 12 


■ o 2 , 










015 




Tl ■ 


T 3 fi 3 


• O i fi3 


■ o 3 , 










016 




T2 ' 


T 3 fi2 


' 02 fu 


■ o 3 , 










017 




Tl ■ 


T 3 f 23 


■ Ol fl2 


■ o 3 , 










018 




Tl ■ 


(T2 X 


T 3 )oi ■ 


o 3 o 2 ■ 


(ri2 x 


f*2 3 ) , 






019 




Tl ■ 


(r 2 x 


T 3 )o 3 ■ 


f 23 f 23 


■ (oi x 


fa) , 






020 




Tl ■ 


(T2 X 


T 3 )oi ■ 


f 23 o 2 


• r 23 o 3 


■ (fl2 


X 


^23) , 


021 




Tl ■ 


(r 2 X 


T 3 ) Ol ■ 


fi 3 o 3 


■ fl 3 fa 


■ (fl2 


X 


^23) , 


022 




Tl ■ 


(T2 X 


T 3 ) O l • 


f 23 o 3 


■ f 12 fa 


■ (fl2 


x 


^23) , 



(6.43) 

where = rij/\rij\ and = fi — fj denotes the position of nucleon i with respect to nucleon j. The 3NF is a linear 
combination of the operators Gi with the coefficients given by scalar functions Ti(r\2, r 23 , r 3 i). These functions have 
the dimension of energy and can be interpreted as the potential energy between three static nucleons projected onto 
the corresponding operator. The profile functions Ti receive contributions from the long-range and the intermediate- 
range 3NF topologies and are predicted (at long distances) in terms of the chiral expansion. In order to explore the 
convergence, we plot these functions for the equilateral triangle configuration of the nucleons given by the condition 

ri 2 = r 23 = r 31 = r . (6.44) 

Restricting ourselves to this particular configuration allows us to stay with simple one-dimensional plots. We em- 
phasize, however, that the conclusions about the convergence of the chiral expansion for the 3NF drawn in this 
section apply to this particular configuration. We begin with the longest-range 2ir exchange topology. Projecting the 
coordinate-space expressions given in section [TT] onto the operators in Eq. (|6.43|) and evaluating the three-dimensional 
integrals in Eqs. (|2.8[) and (|2.13[) numerically we compute the corresponding contributions to the profile functions 
J" (3) (r), J"M(r) and J" (5) (r) at N 2 LO, N 3 LO and N 4 LO, respectively. Our results for the 3NF profile functions 
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FIG. 4: Chiral expansion of the profile functions Ti(r) in MeV generated by the two-pion exchange 3NF topology up to 
N 4 LO (in the equilateral triangle configuration) . Dashed-dotted, dashed and solid lines correspond to , + J-^ and 



•^f ' + ^T' + respectively. 



r (4) 



(5) 



generated by the 2ir exchange diagrams are visualized in Fig. [4] Here and in what follows, we use the values of the 
low-energy constants corresponding to the order-Q 4 KH fit to the pion-nucleon partial wave analysis of our work flol ] . 
In particular, we employ the following values of Ci (all in units of GeV -1 ): 

ci = -0.75, c 2 = 3.49, c 3 = -4.77, c 4 = 3.34. (6.45) 

The results for the functions F%{t) plotted in Fig. U resemble the findings of our work [40| . where a good convergence 
of the chiral expansion in momentum space was observed by looking at the functions A{q2) and $(32) for low values 
of the momentum transfer. While there are large corrections at N 4 LO to some of the profile functions and, especially, 
to Fi(r) at short distances of the order of r ~ 1 fm, we observe a very good convergence at long distances of the 
order of r > 2 fm. At such large distances, the N 4 LO results appear to be very close to N 3 LO ones. As already 
pointed out in the introduction, fast convergence of the longest-range 3NF is not surprising given that effects of the 
A-isobar are, to a large extent, accounted for already in the leading contributions F> (r) to this topology through 
resonance saturation of the LECs C3.4. We further emphasize that the operator structure of the 27r exchange topology 
is fairly restricted: only 10 out of 22 functions ^(r) get non- vanishing contributions. Notice that the larger number of 
nonvanishing functions Fi in coordinate space compared to momentum space has to be expected due to the appearance 
of gradients when performing the Fourier transform. In contrast to the momentum space representation, the number 
of nonvanishing structures in the coordinate space representation of a 3NF is not representative for estimating the 
number of affected nucleon-deuteron polarization observables at a fixed kinematics. 

It is instructive to compare the strength of the three- and two-nucleon potentials. While the long-range three-nucleon 
potentials are considerably weaker than the two-nucleon potentials, they are still not negligible. For example, the 
isovector-tensor and isoscalar central nucleon-nucleon potentials governed by one-pion exchange and (subleading) 
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r [fm] r [fm] 

FIG. 5: Chiral expansion of the profile functions Ti(r) in MeV generated by the two-pion-one-pion exchange 3NF topology up 
to N 4 LO (in the equilateral triangle configuration). Dashed and solid lines correspond to and + J-^ 5 \ respectively. 

two-pion exchange, respectively, have the strength of the order of 3 . . .4 MeV at distances r ~ 2 fm [l5j . These are 
the strongest two-nucleon forces at large distances. The strongest three-nucleon potentials T15 (r) , Tiq (r) and Fn(r) 
reach at such distances the strength of ~ 0.7 ... 1 MeV. We remind the reader that nuclear potentials become scheme 
dependent at short distances below r ~ 1 ... 1.5 fm, where the contributions of short-range topologies start playing 
important role. An estimation of this intrinsic scheme dependence for nucleon-nucleon potentials is provided in Fig. 3 
of Ref. [13. 

The convergence of the chiral expansion for the two-pion-one-pion exchange and ring topologies is, as expected, much 
worse, see Figs. [5] and |6] In nearly all cases, the subleading contributions at N 4 LO dominate over the nominally 
leading ones at N 3 LO even at large distances. This is analogous to the pattern observed for the two-pion exchange 
two-nucleon potential. In that case, the strong dominance of the subleading terms appears because of several reasons 
including the large numerical coefficients, an enhancement by one power of it as compared to the standard chiral power 
counting which is characteristic to the triangle diagrams, see also Ref. [Hj|, as well as the large numerical values of the 
LECs C3.4 from the subleading pion-nucleon effective Lagrangian which are governed by the A isobar. In the case of 
the 3NF 27T-17T exchange and ring potentials the situation is less dramatic. In particular, the enhancement by a power 
of 7r actually affects the leading contributions at N 3 LO which involve the loop function A(q2). Still, the corrections 
at N 4 LO are large which can presumably be attributed to the large numerical values of the LECs Cj. One should, 
however, emphasize that the potentials generated by the 2it-1it exchange and ring diagrams have a considerably 
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FIG. 6: Chiral expansion of the profile functions Ti{r) in MeV generated by the ring 3NF topology up to N 4 LO (in the 
equilateral triangle configuration). Dashed and solid lines correspond to and + respectively. 



shorter range as compared to the 2-7T exchange ones and only reach at most ~ 50 keV at distances of the order of 
r ~ 2 fm. It is, therefore, not clear whether the lack of convergence will have any significant phcnomcnological impact. 
Clearly, to answer this question one needs to carry out complete calculations of few- and many-nucleon observables. 
This work is in progress. Last but not least, we emphasize that especially the ring topology generates a very rich 
structure in the 3NF and gives rise to 20 out of 22 profile functions. 
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It is also instructive to compare the 3NF potentials corresponding to the individual topologies with each other. Such 
a comparison is given in Fig. where we restrict ourselves to N 4 LO, i.e. we only show jf } + jf ) + jf } . We 
observe that the 2ir-lir exchange and ring potentials arc of a comparable size. However, in all cases where the longest- 
range 2ir exchange topology contribute, it clearly dominates at r > 2 fm over the two other topologies. At shorter 
distances of the order of r ~ 1 fm the impact of the 2ir-lir exchange and ring terms becomes more significant with, 
e.g. |J r n i i5(l fm)| ~ 20 MeV to be compared with the strongest 2tt exchange potentials | ^-"15,16,17(1 fni)| ~ 100 MeV. 
As pointed out before, it is difficult to draw conclusions on the phcnomcnological importance of the new structures 
based on this comparison alone since one generally expects that (scheme-dependent) short-range contributions to the 
3NF not considered in the present work would become significant at r < 1 fm. 

Last but not least, Fig. [8] shows the chiral expansion of the complete results for the functions Fi(v) emerging from 
adding the contributions from all three topologies together. The interpretation follows directly from the above 
discussion. At long distances of the order of r > 2 fm dominated by the 2tt exchange one observes a good convergence 
for all cases where the potentials are numerically sizable. On the other hand, those profile functions which are not 
affected by the 2ir exchange are typically dominated by the N 4 LO contributions which might be still not converged 
at this order in the chiral expansion. The corresponding potentials are, however, rather weak. At shorter distances 
j- ~ 1 ... 2 fm, the 2ir-lir exchange and ring contributions start becoming more important relative to the 2tt exchange 
terms. One again observes the dominance of the N 4 LO contributions which supports the assumption about the 
important role played by A excitations, whose effects are partially taken into account at N 4 LO through resonance 
saturation of the LECs C2, C3 and C4. 



VII. SUMMARY AND OUTLOOK 



In this paper, we have worked out and analyzed in detail the intermediate-range contributions to the three-nucleon 
force at N 4 LO, which emerge from the 2-k-Itt exchange and ring topologies. We used here the heavy-baryon formulation 
of chiral EFT with pions and nucleons being the only explicit degrees of freedom. The pertinent results of our study 
can be summarized as follows. 

• We worked out the coordinate-space representation of the N 4 LO corrections to the 2tt exchange 3NF calculated 
in momentum space in Ref. [40| | . 

• We derived the N 4 LO contributions to the intcrmcdiatc-rangc 2ir-lir exchange and ring topologies. These 
represent the first corrections to the leading results which appear at N 3 LO and have been worked out in 
Ref. [2^ |. We provide explicit analytical expressions in both momentum and coordinate spaces. The obtained 
corrections do not involve any unknown low-energy constants. 

• We have demonstrated that the most general structure of an isospin-invariant local three-nucleon force involves 
89 independent isospin-spin-momentum operators. We proposed a set of 22 linearly-independent operators which 
can serve as a basis and gives rise to all 89 structures in the 3NF upon making permutations of nucleon labels. 
We also discussed the properties of the corresponding scalar structure functions J-1...22 parametrizing the 3NF 
with respect to permutations. 

• Finally, using the above mentioned operator basis, we addressed the convergence of the chiral expansion for the 
long-range tail of the 3NF in the equilateral triangle configuration with ri2 = ^23 = r^i = r by comparing our 
predictions for the potentials Ti at different orders in the chiral expansion. Consistently with the momentum- 
space results of Ref. we observe a good convergence for the longest-range 2-7T exchange topology which 
clearly dominates the 3NF at distances of the order r > 2 fm. The intermediate-range 2ir-lir exchange and ring 
diagrams provide sizable corrections to J 7 ; at r ~ 1 fm and also contribute to those 12 profile functions which 
vanish for the 2ir exchange. As expected, we found that N 4 LO corrections to the intermediate-range topologies 
are numerically large and in most cases dominate over the nominally leading N 3 LO terms. This can be traced 
back to the role played by the A(1232) isobar whose excitations provide an important 3NF mechanism. In the 
standard, delta-less formulation of chiral EFT we employ here, effects of the A isobar are not incorporated in 
N 3 LO contributions to the 3NF. For the intcrmcdiatc-rangc topologies we are primarily interested in here, first 
effects of the A appear at N 4 LO through resonance saturation of the LECs C2, C3 and C4 which accompany the 
sublcading pion-nucleon vertices in the effective Lagrangian. The importance of the A isobar is reflected in the 
numerically large values of these LECs which are responsible for large N 4 LO corrections we observe. 
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FIG. 7: Individual contributions of the two-pion exchange (dotted lines), two-pion-one-pion exchange (long-dashed lines) 
and ring (dashed-double-dotted lines) topologies to the profile functions Fi{r) in MeV at N 4 LO in the equilateral triangle 
configuration. 
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equilateral triangle configuration). Dashed-dotted, dashed and solid lines correspond to J-^ , J-f^ 1 + J-"; 4 ' and J-"^ 3 ' + J-^ + J-f^ 
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The results of our work provide important step towards precise, quantitative theoretical description of the 3NF in the 
framework of chiral EFT. The long-range part of the 3NF is governed by exchange of pions, the Goldstone bosons 
of the spontaneously broken chiral symmetry of QCD, and can be rigorously calculated in the framework of chiral 
EFT. It is expected to affect the energy dependence of the nucleon-deuteron scattering amplitude at low energies 
and might be responsible for the long-standing puzzles such as e.g. the A y puzzle in elastic three- and four-nucleon 
scattering Q. Although the resulting intermediate-range potentials are significantly weaker than the 2ir exchange 
terms, the appearance of new structures might lead to large effects in certain nucleon-deuteron scattering obscrvables. 
It would be interesting in the future to explore this possibility in a systematic way. Clearly, the N LO corrections 
to the short-range part of the 3NF should also be worked out. This work is in progress. Notice that subleading 
contributions to the three- nucleon contact interactions at N 4 LO are discussed in Ref. [45j]. Finally, the large N 4 LO 
corrections for the intermediate-range terms raise an obvious question in regard to whether the chiral expansion for 
these quantities can be truncated at this order. One should especially keep in mind that while the obtained N 4 LO 
corrections do include some of the l/(mA — itin) contributions through values of the LECs 02,3,4 and, in this sense, 
take into account physics associated with intermediate excitation of a single A isobar, double and triple A excitations 
start contributing only at orders N 5 LO and N 6 LO, respectively. Phenomenological studies of Ref. |39j indicate that 
at least double A excitations might induce sizable 3NFs. This issue must be investigated in the future. Rather than 
calculating N 5 LO and N 6 LO corrections to the 3NF in the delta-less formulation of chiral EFT, which correspond to 
the two-loop level, it is more feasible and probably also more efficient to include the A isobar as an explicit degree 
of freedom in the effective Lagrangian, see Refs. [35l . l36l I46H48I ] for some promising steps in this directions. In such 
a delta-full formulation, the leading types of 1/(toa — rn^), 1/(ma — mN) 2 , ■ ■ ■ contributions are resumed and the 
3NF mechanisms associated with single, double and triple intermediate A-cxcitations arc taken into account already 
at N 3 LO. Work along these lines is underway 
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Appendix A: Ring contributions in momentum space 



In this appendix we give the lengthy expressions for subleading contributions to ring diagrams in momentum space. 
We employ here the general parametrization of local three-body-force: 

22 

King = y^gj fc»(gi, g3, z), (A.l) 

i=l 

where Gi, • ■ ■ ,^22 are the spin-isospin-momentum operators which we defined in Table [IJ The N LO contributions 
to the structure functions IZi (qi , (73 , z) with z — qi ■ q 3 proportional to g\ are given by 

^ 5) ' 5 ' 4 = T^Jh ^ : °> 935 ^ + C3 > ( 144M * ^ x ) 2 + 8M " ^ - x ) ( 16z2 («? + q ^ 

b\2k2 (z* — 1) v 



512F6 (z 2 - 1) 

-13 (qf + qj) + 16gi<2 3 2 3 - lOgi^-z) + qf (32z 4 - 32z 2 + 9) + Aqlq 3 z (l - 4 
+2q\ql (80z 4 - 78z 2 + 25) + Aq^z (l - 4z 2 ) 2 + qj (32z 4 - 32z 2 + 9)) 



122887r 2 ^6 gi ( z 2 _ x f 



L(q 3 )(c 2 + c 3 ) (4M 2 (z 2 - 1) ( qi (I28z 2 - 101) + 27^) 



-3qf (32z 4 - 32z 2 + 9) + iq\q s z (32z 4 - 5) + qi ql (320z 4 - 412z 2 + 173) + 3g|z (26z 2 - 17)) 

- g j} , 2 r] , L(qi)(c2 + c 3 ) (4Af 2 (z 2 - 1) (27gi2! + q 3 (l28z 2 - 10l)) 
12288^^°^ \z A - 1) 



23 



+g 3 (78z 3 - 51z) + q\q 3 (320z 4 - 412z 2 + 173) + 3q x q$z (32z 4 - 5) + 3gf (32z 4 - 32z 2 + 9)) 

- gg 2w f A , 2 7^ HQ2)(C2 + c 3 ) (-4M 2 (z 2 - 1) (27. (g 2 + g 3 2 ) - 74< M3 * 2 + 128 gi g 3 ) 
122887r 2 F£gig 3 (z 2 - 1) 

+g 4 (51z - 78z 3 ) + iqfqa (2z 4 - 85z 2 + 56) + 2g 2 g 2 z (ll6z 4 - 424z 2 + 227) 
+Aq x q! (2z 4 - 85z 2 + 56) + 3g 4 z (17 - 26z 2 )) - ^ + ^ {q! + ^ + 4) 



4096tt 2 F6 (z 2 - 1) 



^2 5) ^ = ocg j\ - , 2 J(4 : 0, -gi, g3i 0) (32 Cl (z 2 - l) (4 (3A/ 2 + q\ + g 2 ) z 2 + 2 gl g 3 z 
256i< ° (z 2 — 1) 

-3 (4A/ 2 + g 2 + g 2 )) A/ 2 + c 3 (-240 (z 2 - l) 2 A/ 4 - 24 (z 2 - l) (8 gi g 3 z 3 + 8 (g 2 + q 2 ) z 2 

-6gig 3 z - 7 (g 2 + g 2 )) A/ 2 + g 4 (-32z 4 + 48z 2 - 19) + g 4 (-32z 4 + 48z 2 - 19) 
-4g x g 3 z (I6z 4 - 20z 2 + 7) - 4g 3 g 3 z (l6z 4 - 20z 2 + 7) - 2g 2 g 2 (64z 4 - 90z 2 + 35))) 

+ L(qi) 

6144F67T 2 ( 4M 2 + g2) qs {z 2 _ 1} 2 ( _ 4 (z 2 _ !) M 2 + g 2 + + 2(?1(?3Z) 

x (96ci (z 2 - 1) (16 (z 2 - 1) (3q x z + q 3 (8z 2 - 5)) Af 4 + 4 (z (5z 2 - 6) qf + q 3 (l2z 4 - 29z 2 + 14) q\ 
+3g 2 z (3 - 4z 2 ) Ql + gf (3 - 4z 2 )) A/ 2 - g 2 (zg 3 + q 3 (Uz 2 - ll) g 2 + g 2 z (l6z 2 - 13) q x 

+ql (4z 2 - 3))) A/ 2 + c 3 (-64 (z 2 - l) 2 (45 9l z + q 3 (256z 2 - 211)) M% 

-16 (z 2 - 1) (3z (53z 2 - 57) q\ + q 3 (848z 4 - 1673z 2 + 789) q\ 

+4qlz (48z 4 - 152z 2 + 95) q x + 4gf (48z 4 - 124z 2 + 73)) A/ 4 - 4 (3z (46z 4 - 110z 2 + 61) gf 

+q 3 (688z 6 - 2546z 4 + 2782z 2 - 969) q\ + 2q\z (l44z 6 - 992z 4 + 1357z 2 - 554) q\ 

+2gf (48z 6 - 640z 4 + 929z 2 - 382) q\ - Zq\z (96z 4 - 128z 2 + 47) q x 

-3g| (32z 4 - 48z 2 + 19)) A/ 2 + g 2 (g 2 + 2q 3 z qi + q 2 3 ) ((42z 3 - 33z) q\ 

+g 3 (496z 4 - 884z 2 + 415) q\ + 3g 2 z (32z 4 - 32z 2 + 9) qi + 3gf (32z 4 - 48z 2 + 19)))) 

— 5 L(q 2 ) 

6144F67r 2 gig 3 (4A/ 2 + q 2 + q 2 + 2 qi q 3 z) (z 2 - l) 2 (-4 (z 2 - 1) M 2 + q 2 + q 2 + 2q x q 3 z) 

x (c 3 (-64 (z 2 - l) 2 (-166gig 3 z 2 + 45 (g 2 + g 2 ) z + 256gig 3 ) M% + 16 (z 2 - l) (3z (57 - 53z 2 ) g 4 

+4g 3 (53z 4 - 281z 2 + 240) g 3 + 2g 2 z (434z 4 - 1433z 2 + 1035) q\ + 4g| (53z 4 - 281z 2 + 240) gi 
+3g|z (57 - 53z 2 )) M 4 - 4 (g 2 + 2g 3 zgi + g 2 ) (3z (46z 4 - 110z 2 + 6l) g 4 
-4g 3 (34z 6 - 314z 4 + 577z 2 - 288) g 3 + 2g 2 z (-268z 6 + 1574z 4 - 2476z 2 + 1143) q\ 
-4g| (34z 6 - 314z 4 + 577z 2 - 288) q x + 3g^z (46z 4 - 110z 2 + 6l)) M 2 
+ {q\ + 2g 3 zgx + g 2 ) 2 ((42z 3 - 33z) qj - 4g 3 (82z 4 - 203z 2 + 112) q\ 
-2q\z (364z 4 - 824z 2 + 433) q\ - 4gf (82z 4 - 203z 2 + 112) q x + 3q 3 z (Uz 2 - ll))) 
-96ciM 2 (z 2 - 1) (-16 (z 2 - 1) (3zg 2 + 2g 3 (8 - 5z 2 ) q x + 3g 2 z) A/ 4 + 4 (z (6 - 5z 2 ) qj 
+4g 3 (z 4 - 8z 2 + 8) q\ + 2q 2 3 z (lOz 4 - 41z 2 + 34) g 2 + 4g 3 (z 4 - 8z 2 + 8) gi + q\z (6 - 5z 2 )) A/ 2 

+ (g 2 + 2g 3 zgx + g 2 ) 2 [zq\ + 2q 3 (8 - 7z 2 ) q x + g 2 z)) ) 

+ 7 T~~ L{q 3 ) 

6144F^7r 2 gi (4Af 2 + g|) (z 2 - 1) (-4 (z 2 - 1) A/ 2 + g 2 + g| + 2gig 3 z) 

x (96ci (z 2 - 1) (16 (z 2 - l) (3g 3 z + q x (8z 2 - 5)) A/ 4 - 4 ((4z 2 - 3) g 3 + 3g 3 z (4z 2 - 3) g 2 

+g 2 (-12z 4 + 29z 2 - 14) q x + g 3 z (6 - 5z 2 )) M 2 - g 2 ((4z 2 - 3) g 3 + g 3 z (l6z 2 - 13) g 2 

+g 2 (I4z 2 - 11) q x + g 3 z)) Af 2 + c 3 (-64 (z 2 - l) 2 (45g 3 z + q x (256z 2 - 211)) AfJ 

-16 (z 2 - 1) (4 (48z 4 - 124z 2 + 73) g 3 + 4g 3 z (48z 4 - 152z 2 + 95) q\ 



24 



+gf (848z 4 - 1673z 2 + 789) gi + 3gfz (53z 2 - 57)) Af 4 + 4 (3 (32z 4 - 48z 2 + 19) q\ 
+3g 3 z (96z 4 - 128z 2 + 47) qf + 2g| (-48z 6 + 640z 4 - 929z 2 + 382) qf 
+2gfz (-144z 6 + 992z 4 - 1357z 2 + 554) q\ + qj (-688z 6 + 2546z 4 - 2782z 2 + 969) qi 
-3qlz (46z 4 - 110z 2 + 61)) M 2 + q 2 (q\ + 2q 3 z qi + q 2 ) (3 (32z 4 - 48z 2 + 19) q\ 
+3g 3 z (32z 4 - 32z 2 + 9) g 2 + gf (496z 4 - 884z 2 + 415) q 1 + 3gf z (l4z 2 - 11)))) 
c 3 (g 2 + 2q 3 z qi + gf ) g\ 
2048F> 2 (z 2 - 1) ' 

ftf'^ = - 16f6 3 ( g J_ x) ^ 4 : 0.-gi.g3;0)cagig 3 (gii8 + (fe)(gi +<fag) 



3ff 



-giL(gi) (8ciA/ 2 (z 2 - 1) (g! + g 3 z) 



128tt 2 F6 (z 2 - 1) (-4M 2 (z 2 - 1) + q\ + 2q lQ3 z + q 2 ) ' 

+c 3 (-8A/ 2 (z 2 - 1) (g! + g 3 z 3 ) - (q\ + 2q iq3 z + gf) (2q x z 2 - 3gi - ga*))) 

3g\ 

+ 128^ ( Z 2 _ 1} (_ 4M 2 (z 2 -l) + q l + 2qi q 3 z + q 2 ) ' 

+c 3 ((g 2 + 2gig 3 z + g 3 2 ) (giz - 2<? 3 z 2 + 3g 3 ) - 8A/ 2 (z 2 - l) ( Ql z 3 + g 3 ))) 

3fll 

128tt 2 F6 (z 2 - 1) (4M 2 + g 2 + 2g!g 3 z + gf ) (-4M 2 (z 2 - 1) + q\ + 2q lQ3 z + q 2 ) 

X (8 Cl A/ 2 (z 2 - 1) (4Af 2 (g x ( Ql - q 3 z 3 + 3g 3 z) + q 2 3 ) + {q\ + 2g x g 3 Z + gf) 2 ) 

+C 3 (-32M 4 (z 2 - 1) (ql - q m z (z 2 - 3) + g 2 ) - 4M 2 (q 2 + 2 gi q 3 z + q 2 ) (q 2 (4z 2 - 5) 

+gig 3 z (-3z 4 + 15z 2 - 14) + ql (Az 2 - 5)) - (g 2 + 2g!g 3 Z + q 2 ) 2 (2z 2 (ql + q 2 3 ) - 3 (qj + g 2 ) 

+5gig 3 z 3 - 7gig 3 z))) , 



g 3 L(g 3 ) (8ciM 2 (z 2 - 1) (q lZ + q 3 ) 



-L(ga) 



T^'"* = 8f 6 _ n J ( 4 : °> ~gl. 0)glg 3 (c 2 + C 3 )(giZ + g 3 )( gl + q 3 z) 

9a rt \i , \( 2 ( 2 o\ , 2\ 5a?i L (9i)( c 2 +c 3 )(gi + q 3 z) 

L{q2){C2 +c 3 ) (gf -q 1 q 3 z(z -3j+g 3 j 



64tt 2 F6 (z 2 - 1) wyv z ' J/ vyi v ' ' 64tt 2 F6 (z 2 - 1) 

3lg3^(g3)(c2 + c 3 )(giz + q 3 ) 



64tt 2 F6 (z 2 - 1) 



4 

^5 5) ' 5i = - 8F6 ^ _ jr J(4 : 0,-gi,g 3 ;0)c 4 giz( gl z + g 3 ) (q 2 + 2 gi q 3 z + q 2 ) 



„4 

C 4 gii(gi) ((g 2 + 2g ig3 z + g 2 ) (z ( g 2 + gf) 



64tt 2 F6 (?3 (z 2 - 1) (-4M 2 (z 2 - 1) + g 2 + 2 9l93 z + g 2 ) 
+4 (M3 z 2 - 2^x93) - 8A/ 2 (z 2 - 1) (<z 2 z + q lQ3 (3z 2 - l) + g 2 z 3 )) 



c 4 i(g 2 ) (8M 2 (z 2 - 1) ( g 3 z + qjq 3 (z 2 + 2) 



647r 2 F6 g3 (z 2 - 1) (-4M 2 (z 2 - 1) + g 2 + 2 g i g3 z + gf) 
-gig 2 z (z 2 - 4) + g 3 ) - (ql + 2g 1(?3 z + g 2 ) 2 (g x z + g 3 (3 - 2z 

ffl 4 2 

+ 64tt 2 F6 (z 2 - 1) (4M 2 + q 2 ) (-4M 2 (z 2 - 1) + q 2 + 2 gi g 3 z + q 2 )^^ (" 32M - ~ ^ 
X ( 9 2 z 2 + gi93 z (z 2 + 1) + gf) + 4M 2 (q\z 2 + q\q 3 (3z 3 + z) + q\q 2 (-Z 4 + 5z 2 + 2) 
-g!g 3 z (3z 4 + z 2 - 8) + q 4 3 (5 - 4z 2 )) + gf ( 9 2 + 2q iq3 z + gf) (g 2 z 2 + qi q 3 z (z 2 + l) 
+g 3 2 (3-2z 2 ))), 

n ^^ A = oor , fi ^ ; 3^(4 : 0, -gi, g 3 ; 0)c 4 z (q 2 + 2q 3 z gi + q 2 3 ) (4 (3M 2 + q 2 + gf ) z 2 + 2g 1(73 z 

128F°gi(? 3 (z 2 - 1) 



-3(4Ml+q 2 +q 2 )) 



3072f*7T 2 «i (4M 2 + q 2 ) q 2 (z 2 - if (-4 (z 2 - 1) M 2 + q\ + q 2 + 2 qi q 3 z) 
x ( - 1024( 9 i + q 3 z) (z 2 - l) 2 M® - 16 (z 2 - 1) ((41z 2 - 48) q\ + q 3 z (50z 2 - 71) q\ 
+g 3 (24z 4 - 29z 2 - 16) qx + q\z (2Az 2 - 31)) M* - 4 ((25z 4 - 70z 2 + 42) g? 
+g 3 z (40z 4 - 143z 2 + 88) qf + q\ (36z 6 - 115z 4 + 17z 2 + 32) q\ + q\z (l2z 4 - 101z 2 + 59) q\ 
+3g 4 z 2 (7 - 12z 2 ) qi + 3g|z (3 - 4z 2 )) M 2 + q\ (q\ + 2g 3 z 9l + q 2 3 ) ((l3z 2 - 10) q\ 
+q 3 z (I6z 2 - 7) q\ + Zqlz 2 (Az 2 - l) g x + 3gfz (4z 2 - 3))) 

4 

■c A L{q 2 ) (256 (g 2 + 2g 3 zg! + qi) 



mi2F^ 2 q\q 2 (z 2 - if (-4 (z 2 - 1) Ml + g 2 + g 2 + 2 qi q 3 z) 
x (z 2 - l) 2 A/ 4 + 4 (z 2 - l) ((25z 2 - 32) g 4 + 4g 3 z (l3z 2 - 20) q\ + 2gf (l4z 4 - 15z 2 - 20) g 2 
+4gfz (I3z 2 - 20) qx + g 4 (25z 2 - 32)) M 2 - (qj + 2q 3 zq x + q 2 f {2 qi q 3 z 3 + 13 (g 2 + g 2 ) z 2 
+4 qi q 3 z- 10 (q 2 +ql))) 



-CiL{q 3 ) 



Ql 



o 1 

3072FS^q 3 (4Af 2 + g 2 ) (z 2 - l) 2 (-4 (z 2 - 1) M 2 + q 2 + q 2 + 2 qi q 3 z) 

x ( - 1024(g 3 + qxz) (z 2 - l) 2 M% - 16 (z 2 - l) (z (24z 2 - 3l) g? + g 3 (24z 4 - 29z 2 - 16) q\ 

+q 2 3 z (50z 2 - 71) q x + gf (41z 2 - 48)) M 4 + 4 (3z (4z 2 - 3) qf + 3g 3 z 2 (l2z 2 - 7) g 4 

+q\z (-12z 4 + 101z 2 - 59) q\ - gf (36z 6 - 115z 4 + 17z 2 + 32) q\ + g 4 z (-40z 4 + 143z 2 - 8 

+gf (-25z 4 + 70z 2 - 42)) A/ 2 + g 2 (g 2 + 2g 3 zgx + g 2 ) (3z (4z 2 - 3) q\ + 3g 3 z 2 (4z 2 - l) g 2 

_l 2 / 1R 2 7 \ , j/„ 2 i n\ \ \ c 4 (2 gi g 3 + (g 2 + gf) z) ff 4 
+g 3 z(!6z -7)g 1+ g 3 (l3z - 10))) - -^———— r , 

—-, ^ -2 ^(4 : 0, -g l5 g 3 ; 0)g 3 (c 3 (4z 3 (3A'/ 2 + g 2 + g 2 ) - 12Af 2 z - 7z (g 2 + g 3 2 ) + Sg^z 4 
16F°gi (z l - 1) 

-12g!g 3 z 2 - 2g ig3 ) - 16ciA/ 2 z (z 2 - l)) 

3 — L(qx) (cs (4A/ 2 (z 2 - 1) ( 9l (z 2 + 4) 

1287r 2 ^V(z 2 -l) 2 (-4A4 2 (z 2 -l) + g 2 + 2g 1 g 3 * + g§) WA V " V n " 1V j 

+5g 3 z) + (g 2 + 2 gi g 3 z + g 2 ) (3g x (z 2 - 2) + g 3 z (4z 2 - 7))) - 16 Cl A/ 2 (z 2 - l) (qx + g 3 z)) 

q 3 L(q 3 ) (c S {AM 2 (z 2 - 1) (5qxz 



1287r 2 F6g 2 (z 2 - l) 2 (-4M2 (z 2 - 1) + g 2 + 2gxg 3 z + g 2 ) 
+g 3 (z 2 + 4)) + (g 2 + 2gxg 3 z + g 2 ) (q lZ (Az 2 - 7) + 3g 3 (z 2 - 2))) - 16 Cl A/ 2 (z 2 - l) ( 9l z + g 3 )) 

5 — L(q 2 ) 

1287T 2 F^q 2 (z 2 - if (AMI + q\ + 2<7i^ + q 2 3 ) (-4A/ 2 (z 2 - 1) + g 2 + 2g ig3 z + g 2 ) 

x (l6 Cl Af 2 (z 2 - 1) (4Af 2 (qx (qx - g 3 z 3 + 3g 3 z) + g 2 ) + (g 2 + 2g x g 3 z + g 2 ) 2 ) 

+c 3 (-16M 4 (z 2 - 1) (g 2 (z 2 + 4) - 2gxg 3 z (z 2 - 6) + g 2 (z 2 + 4)) - 4A/ 2 (g 2 + 2gxg 3 z + g 2 ) 
x (g 2 (z 4 + 6z 2 - 10) - 2gig 3 z (2z 4 - 13z 2 + 14) + q\ (z 4 + 6z 2 - 10)) 

- (q\ + iqiqsz + <&f (Ml (z 2 - 2) + 2g l93 z (4z 2 - 7) + 3g 2 (z 2 - 2)))) + 

2 /(4 : 0, -gi, g 3 ; 0) (l6 Cl A// 2 (z 2 - l) + c 3 (-12A/ 2 (z 2 - l) - 2z 2 (g 2 + 3 qi g 3 z + g 2 ) 

32-^^ (z 1) 

+5 (g 2 +g 2 ) +12gxg 3 z)) 
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256w 2 F%q 3 (z 2 - if (—AM 2 (z 2 - 1) + g 2 + 2gxg 3 z + g 2 ) 
+c 3 (-4A/ 2 (z 2 - 1) (5g lZ + g 3 (4z 4 - 4z 2 + 5)) - (g 2 + 2g x g 3 z + g 2 ) 



L(q x ) (WcxM 2 z (z 2 - l) ( Ql + g 3 z) 
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x {A qi z 3 - 7q lZ + 2q 3 z 2 - 5q 3 ))) 
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2567r 2 F6 gi {z 2 - \f {-AMI ( z2 - I) + q 2 + 2q x q 3 z + q 2 ) 



L(q 3 ) {l6 Cl M 2 z (z 2 - 1) {q x z + q 3 ) 



+c 3 {-AMI (z 2 ~ 1) (5i (4^ 4 - 4z 2 + 5) + 5g 3 z) - (<ff + 2<7ig 3 z + g§) 
x (2 9l z 2 - 5gi + 4g 3 z 3 - 7g 3 z))) 
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256n 2 F% qi q 3 (z 2 - l) 2 (4A/ 2 + q\ + 2q lQ3 z + q 2 ) {-AMI (z 2 - 1) + qf + 2 Ql q 3 z + q 2 ) 
x (16 C1 M 2 (z 2 - 1) (AMI (qi (qi - q 3 z 3 + 3q 3 z) + q 2 ) + {q 2 + 2 gi q 3 z + q 2 ) 2 ) 
+C 3 (-16M* (z 2 - 1) (5 (q 2 + q 2 ) - Aq x q 3 z 3 + Uq iq3 z) - AM 2 (3z 2 - 4) {q 2 + 2q lQ3 z + q 2 ) 
X (3 {q\ + q^ - 2q lQ3 z 3 + 8 qi q 3 z) - {q\ + 2 Ql q 3 z + q 2 ) 2 {Az 2 {q 2 + q 2 ) - 7 {q\ + qj) 

+10q iq3 Z 3 I6q iq3 z))) ^ 



256n 2 F%{z 2 - 1)' 



n ^' A = aoJl 9 * iv^ 4 : °' ~ qi ' qs; ° )Z ( 16c ^ M " z ( z2 ~ !) + c 3 H* 3 ( 3M ' + + «§) + 12M ' Z 
32-Ttt (z 1) 

+7z (<? 2 + g 2 ) - 8<?i<? 3 z 4 + \2 qi q 3 z 2 + 2q 1 q 3 )) 



2567r 2 F6g 3 (z 2 - l) 2 (-4M 2 (z 2 - 1) + g 2 + 2?ig&* + g 2 ) 
+c 3 (-4M 2 (z 2 - 1) { Ql (z 2 + 4) + 5g 3 z) - {q 2 + 2q x q 3 z + q\) (3<ft (z 2 - 2) + g 3 z (4z 2 - 7)))) 



zL(gi) (16 C1 M 2 (z 2 - 1) (<zi+g 3 z) 
g 3 2 ) (3gi (z 2 -2)+g 3 z(4z 2 -7))); 
zL{q 3 ) (16 C1 M 2 (z 2 - 1) {q lZ + q 3 ) 



256ir 2 F%q x { z 2 _ i)^ (_ 4M 2 (z 2 _ l) + 9 2 + 2<M3 z + g 2 ) 
+c 3 (-4A/ 2 (z 2 - 1) {5q x z + q 3 (z 2 + 4)) - {q 2 + 2 Ql q 3 z + q 2 3 ) { Ql z {Az 2 - 7) + 3q 3 (z 2 - 2)))) 



2567T 2 F6g ig3 (z 2 - 1)" (4Af 2 + q\ + 2q iq3 z + q 2 ) {—AM 2 (z 2 - 1) + q\ + 2q x q 3 Z + q 2 ) 
x (16ciM 2 (z 2 - 1) (4Af 2 { Ql ( qi - q 3 z 3 + 3q 3 z) + g 2 ) + {q\ + 2q x q 3 z + q 2 ) 2 ) 

+c 3 (-16M* (z 2 - 1) (g 2 (z 2 + 4) - 2q ig3 z (z 2 - 6) + g 2 (z 2 + 4)) - 4A/ 2 (<j 2 + 2q x q 3 z + 
x (g 2 (z 4 + 6z 2 - 10) - 2<7iq 3 z (2z 4 - 13z 2 + 14) + q\ (z 4 + 6z 2 - 10)) 

- {q 2 + 2 gi q 3 z + q 2 ) 2 {iq 2 (z 2 - 2) + 2 gi q 3 z {Az 2 -7)+ 3q 2 (z 2 - 2)) 



256tt 2 F6 (z 2 - 1) ' 



^n'^ = - — K gA 9 ^(4 : 0, -qi,q 3 ; 0)c 4 z{q lZ + q 3 ) {-12M 2 (z 2 - l) + iq\ + 6 gi q 3 z + ql (7 - 4z 2 )) 

lol<°qi {z z — 1) 

51 c 4 i(, 1 )(4Af 2 (z 2 -l)(< ? 2 z(z 2 +4) 



1287^^3 (z 2 - l) 2 (-4A/ 2 (z 2 - 1) + q\ + 2 gi q 3 z + q 2 ) 
-2q lQ3 (2z 4 - 9z 2 + 2) + 5q 2 z) - {q 2 + 2q lQ3 z + q 2 ) {q\z (z 2 + 2) - 2 gi93 (2z 4 - 7z 2 + 2 
+g 3 2 z(7-4z 2 ))) 



Ci L{q 2 ){{q 2 + 2q iq3 z + q 2 ) 



128n 2 F%q 2 q 3 (z 2 - l) 2 (-4A/ 2 (z 2 - 1) + q 2 + 2 Ql q 3 z + q 2 ) 
x {q\z (z 2 + 2) + 3<? 2 g 3 (2z 4 - z 2 + 2) - q x q\z (z 2 - 10) - 3gf (z 2 - 2)) 
-4M 2 (z 2 - 1) ( z 2 + 4) + q\q 3 (6z 4 + z 2 + 8) + 3 M 2 z (z 2 + 4) + q\ (z 2 + 4))) 

5 ^ CiL{q 3 ) (16A/ 4 (z 2 - 1) 

^Stt 2 ^ 2 ( z 2 _ !) 2 ( 4M 2 + g |) ( _ 4M 2 ( z 2 _ 1) + g 2 + 2qi q 3 Z + q 2 ) M 1 7 

x (5g 2 z 2 + 5qi(7 3 z (z 2 + l) + g 2 (z 2 + 4)) - 4A/ 2 (3q 4 z 2 + 3q 3 q 3 (3z 3 + z) + g 2 q 2 (-z 4 + 15z 2 + 4) 
+ Ql qlz (-7z 4 + z 2 + 18) - ql (z 4 + 6z 2 - 10)) - 3q 2 {qj + 2q x q 3 z + qj) {qjz 2 + q x q 3 z (z 2 + l) 
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128^ 2 F6 gi (z 2 -l)' 

n ^' A = ~ 1CFfi ^2 ^2 J ( 4 : °> 93i 0)c 4 2 (l2M 2 (z 2 - 1) ( 9l + q 3 z) + q\ (- (2z 2 + l)) 
lbb%q 3 (z l — 1) 

-q\q 3 Z (4z 2 + 5) - ftft 2 (4z 2 + 5) + g 3 Z (4z 2 - 7)) 

-c 4 zL{q 1 )((q 2 1 +2q 1 q 3 z + ql) 



U8n*FSql (z 2 - i f (-4M| (z 2 - 1) + ft 2 + 2ftft* + ft 2 ) " 
X (3ft 2 z + qi (5ftz 2 + q 3 ) + q 2 z (7 - 4z 2 )) - 4M 2 (z 2 - l) (5z (ft 2 + ft 2 ) + 9ft ftz 2 + qi q 3 )) 

n C4zL(ft) (q? + 2q±q 3 z + q 2 ) 

12^F% qi q 2 (z 2 - l) 2 (-4Af 2 (z 2 - 1) + q\ + 2ftftz + ft 2 ) ^ ^> 

X (4M 2 (z 2 - 1) (5ft z + q 3 (z 2 + 4)) - 3ft 3 z - qjq 3 (7z 2 + 2) + qi q 2 z (2z 2 - ll) + 3g| (z 2 - 2)) 

s — c 4 zL(ft) 

128n 2 F®qiq 3 (z 2 — 1) (4A/ 2 + q 2 ) (— 4Af 2 (z 2 — 1) + q 2 + 2qiq 3 z + g|) 

x ( - 16 A/ 4 (z 2 - 1) (ft 2 (4z 2 + 1) + 2 9l q 3 z (2z 2 + 3) + ft 2 (z 2 + 4)) + 4Af 2 (qf (2z 2 + l) 

+6qfq 3 z (z 2 + l) - 2ft 2 ft 2 (z 4 - 7z 2 - 3) - 2ft q\z (3z 4 + z 2 - 10) - q 4 (z 4 + 6z 2 - 10)) 

+ft 2 (ft 2 + 2ftftz + ft 2 ) (ft 2 (2z 2 + 1) + 2ftftz (z 2 + 2) - 3ft 2 (z 2 - 2))) 



1287r 2 F6ft (z 2 - 1) : 

n ^' A = ~ 1CF6 g A i , 2 /(4:0,-ft,ft;0) C4 (ft^ + g3)(l2M 2 (z 2 -l)-ft 2 (2z 2 + l) 
16F°ft (z 2 - 1) 

2ftft(2z 3 + z)+ft 2 (2z 2 -5)) 



C 4 L(ft)((ft 2 +2ftftZ + ft 2 ) 



l2%-K 2 F%q 2 (z 2 - l) 2 (—AM 2 (z 2 - 1) + q\ + 2ftftZ + ft 2 ) 
x (3q 2 z 2 + 2ft 93 (2z 3 + z) + ft 2 (5 - 2z 2 )) - 4M 2 (z 2 - l) (5ftV + 2ft ft (4z 3 + z) 
+ft 2 (4z 4 -4z 2 + 5))) 

3 — c 4 L(q 2 ) UM 2 (z 2 - 1) (5ftV 

1287T 2 F6ftg 3 2 (z 2 - l) 2 (-4M 2 (z 2 - 1) + q\ + 2q iq3 z + ft 2 ) V M 

+5g 2 ft (2z 3 + z) + ftft 2 (-4z 4 + 23z 2 - 4) + 5ft]z) - (ft 2 + 2ftftz + ft 2 ) (3ftV + 3ft 2 © (2z 3 + z) 

+ftg 2 (-8z 4 + 21z 2 - 4) + qjz (7 - 4z 2 ))) 

3 — CiL{q 3 ) (-16Af 4 (z 2 - 1) 

1287r 2 F6ft93^ 2 -l) 2 (4M 2 +ft 2 )(-4Af 2 (z 2 -l)+g 2 + 2ftftz + g 2 ) n V ' 

x (ft 2 (4z 3 + z) + ftft (4z 4 + 5z 2 + l) + 5g 2 z) + 4M 2 (qf (2z 3 + z) + q\q 3 (6z 4 + 5z 2 + l) 

+q\q\z (-2z 4 + llz 2 + 9) + ftft? (-6z 6 - 3z 4 + 19z 2 + 2) + iq\z (4 - 3z 2 )) 

+q 2 3 (ft 2 + 2ftftz + ft 2 ) (ft 2 (2z 3 + z) + ftft (2z 4 + 3z 2 + 1) + q 2 z (7 - 4z 2 ))) 

fflc 4 (ftz + ftO 
1287r 2 F6g 3 (z 2 -l)' 

4 

^14' 944 = - TT^CTh ^2 J ( 4 : 0,-ft,ft;0)c 4 ft (-12M 2 (z 2 - 1) (ft +q 3 z) + qf (2z 2 + 1) 
16F6g 2 (z 2 - 1) 

+g 2 ftz (4z 2 + 5) + q iq 2 3 (4z 2 + 5) + qlz (7 - 4z 2 )) 



-c 4 fti(ft) ((ft 2 + 2ftftz + ft 2 ) 



1287r 2 ^ g 3 (z 2 _ 1} 2 ( _ 4M 2 (z 2 _ 1) + g 2 + 2 ftftZ + ft 2 ) ' 

x (3ft 2 z + ft (5ftz 2 + ft) + q 2 z (7 - 4z 2 )) - 4M 2 (z 2 - l) (5z (ft 2 + ft 2 ) + 9ftftz 2 + ftft)) 



A 



-c 4 L(q 2 ) (q 2 + 2ftftz + qf) 



l2%n 2 F%ql (z 2 - l) 2 (-4A/ 2 (z 2 - 1) + q 2 + 2ftftz + ft 2 ) 
x (-4A/ 2 (z 2 - 1) (5ftz + ft (z 2 + 4)) + 3qfz + ft 2 ft (7z 2 + 2) + q iq 2 z (ll - 2z 2 ) - 3ft 3 (z 2 - 2)) 
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5 — CiL(q 3 ) (16M* (z 2 - 1) 

1287^2 (z 2 _ 1)2 (4M 2 + g 2) ( _ 4M 2 (z 2 _ 1} + q 2 + 3^ + g 2) M V 7 

x (g 2 (4z 2 + 1) + 2 qi q 3 z (2z 2 + 3) + g 2 (z 2 + 4)) - 4A/ 2 (g? (2z 2 + l) + 6^q 3Z ( z 2 + l) 
-2g 2 g 2 (z 4 - 7z 2 - 3) - 2g!gf z (3z 4 + z 2 - 10) - g 4 (z 4 + 6z 2 - 10)) - g 2 (g 2 + 2 gi q 3 z + q 2 ) 



x (g 2 (2z 2 + 1) + 2 qi q 3 z (z 2 + 2) - 3q| (z 2 - 2))) + 



1287r 2 F6 g 2 (z 2 - 1) ' 



nl & Ul = - ,a^j\ i n2 J ( 4 : °. ~gi» g3i 0)(C2 + c 3 ) (2z 2 (-2M 2 + g 2 + g 2 ) + 4M 2 + g 2 

+2g!g 3 z 3 + 4gig 3 z + q 2 3 ) 
g\zL{q 2 ){c 2 + c 3 ) (3g 2 - 2 gl g 3 z (z 2 - 4) + 3g|) g 4 1 L(gi)(c 2 + c 3 ) (3giz + 2g 3 z 2 + g 3 ) 
128^ 2 F6g 1 g 3 (z 2 - l) 2 1287r2F6g 3 (z 2 - l) 2 

+ 5ii(g 3 )(c 2 +c 3 )(2giz 2 + g 1 +3g 3 z) f g\{c 2 +c 3 ) 



128^ 2 F6g! (z 2 - l) 2 128^ 2 F6 (z 2 - 1) ' 

4 

^le'' 1 = ^f— — I /(4:0,-g 1 ,g 3 ;0)z( C2 + C3 )(g 2 + 2g 1 g 3 z + g 2 ) 

8F°gig 3 (z 2 - 1) 

x (-4A/ 2 (z 2 - 1) + 3g!(g! + 2g 3 z) + 2g 2 z 2 + g 2 ) 
g 4 A L{q 2 )(c 2 + c 3 ) (g 2 + 2g lg3 z + g 2 ) (g 2 (z 2 + 2) + 2 gl g 3 z (z 2 + 2) + 3g 2 z 2 ) 

647r 2 F6g2g2 ( z 2 _ !)2 

| g\L{ qi ){c 2 + c 3 ) (gg (z 2 + 2) + g 2 g 3 z (4z 2 + 5) + gl g 2 z 2 (2z 2 + 7) + qjz (2z 2 + l)) 

647r 2 ^g 1 g| (z 2 - l) 2 
3g 4 1 zL(g 3 )(c 2 + c 3 )(gi + g 3 z) (g 2 + 2g!g 3 z + g 2 ) g^(c 2 + c 3 ) (z (g 2 + gj) + 2gig 3 ) 
64^ 2 F6g 2 g 3 (z 2 - l) 2 64n 2 F* qi q 3 (z 2 - 1) 

^i7 hA = ~ 16F6 ^ 2 _ 1} ^(4 : 0, -gi, g 3 ; 0)(ca + c 3 ) (-4M 2 (z 2 - l) + g 2 + 2g l93 z + g 2 ) 

g\zL{q 2 ){c 2 + c 3 ) (g 2 + 2gig 3 z + gg) g 4 L(gi)(c 2 +c 3 )(giz + g 3 ) 
1287r 2 F6 gi g 3 (z 2 - 1) 128^ 2 F6g 3 (z 2 - 1) 

gj£(g3)(c 2 + c 3 )(gi + g 3 z) 

1287^31 (z 2 - 1) 

^(5),^ = n (5),g* = n {$),g\ = ^(5), 9 1 = ^(5),^ = ^(5),^ = Q (A 2) 

In the above expressions, gi and g 3 are always to be understood as the magnitudes of the corresponding three-momenta 
(except in the arguments of the function I), gi = |gi |, g 3 = |g 3 |. Further, the function I(d : Pi,P2,P3',Pi) refers to 
the scalar loop integral 

, , 1 f d d l 1 1 1 1 

( :puP 2 ,p 3 ;P4) - .J ( 2n y( l+p j2_ M 2 + ie ( l+p2 }2_ M 2 +ie Q + p3 y_ M 2 + iev .( l + p ^ +ie -( •) 

This expression involves four- momenta pi . For the case pf = which we are interested in, it can be expressed in terms 
of the three-point function in Euclidean space J (d : pi,p 2 ,p 3 ) 

f d d l 1 1 1 

J (d : Pi,P2,Pz) = / , n w — • (A. 4) 

j (27r)d {l+ ^ )2 + (l + + M2 {l+ ^ )2 + M2 

In particular, the function I (4 : 0, — gi, g 3 ; 0) which enters the expressions for Ri can be written as 

I (4 : 0, -g l7 g 3 ; 0) = ^ J (s : 0, -gi, g 3 ) . (A.5) 



The N LO contributions to the structure functions proportional to g\ vanish for 72-1,3,4,6,8,9,10,15,. ..,22- The nonvan- 
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ishing contributions have the form 

U * Wa = ~ ^jh 7T2 J (4 : 0, -gx, g 3 ; 0) (32d (z 2 ~ l) (4 (z 2 - l) M 2 - g 2 - g 3 2 + 2 (g 2 + g 2 ) - 2 
128-f ° (z z — 1) 

+2q iq3 z) Ml - 8c 3 (2M 2 + q\ + qf + 2 gi q 3 z) (z 2 - l) (-4M 2 - g 2 - g 2 + 2 (2M 2 + q\ + q 2 3 ) z 2 
+2q iq3 z) - c 2 (-4A/ 2 - q\ - q% + 4 (M 2 + g 2 + g 2 ) z 2 + 6gig 3 z) 
X (4 (z 2 - 1) Ml - g 2 - q x (q x + 2q 3 z))) 

L(q x ) (-96ci (z 2 - 1) (giz + g 3 (2z 2 - l)) M 2 + c 2 (-3 (2z 3 + z) g? 



3072F> 2 g 3 (z 2 - l) 2 



+q 3 (8 (z 2 - 5) z 2 + 5) qj - 3g 2 z (4z 2 + 5) q x + 3gf (l - 4z 2 ) + 4Af 2 (z 2 - l) (3g x * + q 3 (8z 2 - 5))) 
+8c 3 (z 2 - 1) (3zg 3 + q 3 (I6z 2 - 7) ql + 3<?fz (2z 2 + l) q x + 3g| (2z 2 - l) 
+Ml(6q lZ + q 3 (28z 2 -22)))) 

L(g 2 ) (96 Cl (z 2 - 1) (zg 2 - 2q 3 (z 2 - 2) q x + q 2 3 z) M 2 



3072^> 2 9l g 3 (z 2 - l) 2 



-24c 3 (2i\/ 2 + ql + q 2 3 + 2 qi q 3 z) (z 2 - l) (zg 2 - 2g 3 (z 2 - 2) q x + q 2 3 z) 
+c 2 ((g 2 + 2q 3 z qi + q 2 ) (4 gi q 3 z 4 + 6 (g? + qj) z 3 - 2q x q 3 z 2 + 3 (ql + qf) z + 16 qi q 3 ) 
-AMI {f - 1) (-I0qiq 3 z 2 + 3 (g 2 + q 2 ) z + 16 ? i<fe))) 

L(q 3 ) (-96 Cl (z 2 - 1) (q 3 z + Ql (2z 2 - 1)) Ml + c 2 ((3 - 12z 2 ) ql 



3072F67r 2 gi (z 2 - l) 2 



-3g 3 z (4z 2 + 5) ql + q 2 3 (8 (z 2 - 5) z 2 + 5) g x - 3gf z (2z 2 + l) + 4A/ 2 (z 2 - l) 

x (3q 3 z + q x (8z 2 - 5))) + 8c 3 (z 2 - l) ((6z 2 - 3) q\ + 3g 3 (2z 3 + z) ql + q 2 (l6z 2 - 7) q x 



, C2 {ql + 2q 3 z qi + q 2 ) g 2 A 



+3g 3 z + Ml (6q 3 z + q x (28z 2 - 22)))) + 



1024F6tt 2 (z 2 - 1) 
(ql + 2q x q 3 z + qj) 

g 2 A c A L(q 2 )(q x z + q 3 ) (ql + 2gig 3 z + qf) _ g A C4q 1 zL(q 1 ) (gf + 2gig 3 z + gf ) 

647r 2 F6g 3 (z 2 - 1) 647r 2 i^g 3 (z 2 - 1) 

g 2 A c A L(q 3 ) (qlz 2 + q x q 3 z (z 2 + l) + qf) 



K-f' 91 = — 6 ( g j - J(4 : 0, -g 1; g 3 ; 0)e 4gi z( gl z + g 3 ) (g 2 + 2 gi g 3 z + g 2 ) 
+ 



64w 2 F% (z 2 _ !) 

2 2 

?4 5),9A - - 7777^ g 1 9 T^ 4 : °' -91,93; 0)c 4 z (g 2 + 2 9l g 3 z + g 2 ) (4A/ 2 (z 2 - l) + g 2 (2z 2 - l) 

64F°gig 3 (z 2 - 1) 



+2g 1 g 3 z + g 2 (2z 2 -l)) 

A 

15367T 2 F6g 2 g 2 (z 2 - l) 2 
+2 gi g 3 z (z 2 + 2) + g 2 (7z 2 - 4)) 



c 4 L(q 2 ) (ql + 2q x q 3 z + q 2 ) (16M 2 (z 2 - l) + g 2 (7z 2 - 4) 



c 4 £(gi) (16A/ 2 (z 2 - 1) (gj + g 3 z) + g 3 (7z 2 - 4) + g 2 g 3 z (lOz 2 - l) 



15367^^! (z 2 - l) 2 
+3gigfz 2 (2z 2 + l) + 3g|z (2z 2 - l)) 



2 

c 4 L(g 3 ) (16M 2 (z 2 - 1) (q x z + g 3 ) + 3g 3 z (2z 2 - l) + 3g 2 g 3 z 2 (2z 2 + l) 



15367r 2 F6g 2 g 3 (z 2 - 1) 
+ gi g 3 2 z (lOz 2 - 1) + g 3 (7z 2 - 4)) + 



g\c 4 (z (ql + qj) + 2q x q 3 ) 



15367T 2 f 6 gig3 ( 2 2 _ ]_) 



2 

^ii ),9i = 1CFfi % 772 7 ( 4 : 0. -91= 935 0) C4 z(g!Z + g 3 ) (3 (g 2 + 2g x g 3 z + g 2 ) - 4M 2 (z 2 - l)) 
16F°gi (z 2 - 1) 
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g\c 4 (z 2 + 2) L{q 2 )(q lZ + g 3 ) (g 2 + 2 gl g 3 z + gj) g 2 A c 4 zL( qi ) (q 2 (z 2 + 2) + 6gig 3 z + gg|) 
128^F6g 2 g 3 (z 2 - if ' 1287^^3 (z 2 - if 

g 2 A c 4 L(g 3 ) (3g 2 z 2 + 3g lg3 z (z 2 + l) + g 2 (z 2 + 2)) gjc 4 (gi + g 3 z) 

128n 2 F%q 2 (z 2 - l) 2 1287r 2 F| gi ( z 2 _ ^ > 

- — — 9 f 9 ^(4 : 0, -gi, g 3 ; 0)c 4 z ((g 2 + 2 gl g 3 z + g 2 ) (2g x z 2 + q x + 3g 3 z) 

16F°g 3 [z* - 1) 

,,.2/2 lW . u 9A^zL{q 2 ) {gj + 2g lg3 z + g 2 ) (3giz + g 3 (z 2 + 2)) 

-4M ff z -1 (gi+g 3 z) — — — 2 

128TT 2 F%q iq 2 (z 2 - 1) 

g 2 A c 4 zL{q x ) (3z (g 2 + g 2 ) + 5gig 3 z 2 + qi q 3 ) 

128ir 2 F%q 2 (z 2 - l) 2 

g\c 4 zL(q 3 ) (g 2 (2z 2 + l) + 2 gi g 3 z (z 2 + 2) + g 2 (z 2 + 2)) ffic 4 (g 1 2 + g 3 ) 



mir 2 F% qi q 3 (z 2 - if 128n 2 F%q 3 (z 2 - 1) ' 

2 

^i^ 9 " = - 1CFfi ^2 7 ( 4 : °> "Si. is; 0)c 4 (giz + g 3 ) ((2z 2 + l) (g 2 + 2g ig3 z + g 2 ) - AM 2 (z 2 - l)) 

16F°g 3 (z 2 - 1) 

3g 2 1 c 4 zL(g 2 )(giz + g 3 ) (g 2 + 2gig 3 z + g 2 ) ffjc^gi) (3g 2 z 2 + 2gig 3 (2z 3 + z) + g 2 (2z 2 + l)) 

1287r 2 F6 gi g| (z 2 - if 1287T 2 F6g 2 (z 2 - l) 2 

g 2 A c 4 L(q 3 ) (g 2 (2z 3 + z) + g ig3 (2z 4 + 3z 2 + l) + 3g 2 z) g 2 c 4 (giz + g 3 ) 



128^^5! g 3 (z 2 - if 128^ 2 F6g 3 ( Z 2 _ i) ' 



,2 



n< $' 9A = ?\ -2/(4:0,-g 1 ,g 3 ;0)c 4 g 1 ((g 2 +2g 1 g 3 z + g 2 ) (2g 1 z 2 + g 1 +3g 3 z) 

16F£g 2 (z 2 - 1) 

128ir 2 F%q$ (z 2 - 1) 
g 2 A c 4 L{q 3 ) (qf (2z 2 + l) + 2 gl g 3 z (z 2 + 2) + g 2 (z 2 + 2)) 

1287T 2 F6g§ ( z 2 _ !)2 

5 2 1 c 4 giL(gi) (3z (g 2 + g|) + 5gig 3 z 2 + gig 3 ) g\c 4 qi{qi + q 3 z) 



128^ 2 F6gf (z 2 - l) 2 1287r 2 F6g2 (z 2 - 1) ' 



(A.6) 



Finally, the N 4 LO contributions to the structure functions proportional to g A vanish for 72-i,3,...,6,8,...,22- The nonva- 
nishing contributions have the form 

n ^' 9 ° A = -^zh ^ /(4 : '-9i,9 3 ;0) (4M 2 (z 2 - 1) - qi { qi + 2q 3 z) - q 2 ) (32 Cl Af 2 (z 2 - l) 
2b%r ° (z z — 1) 

+3c 2 (-4M 2 (z 2 - 1) + qj + 2q ig3 z + g 2 ) - 8c 3 (z 2 - l) (2A/ 2 + g 2 + 2 gi g 3 z + g 2 )) 

~ g 2F6 1 71 n 2 ^(g2) (<7? + 2grg 3 z + g 2 ) (-32 Cl A/ 2 (z 2 - l) + c 2 (20Af 2 (z 2 - l) 
20487HF°gig 3 (z 2 - 1) 

+ (2z 2 - 5) (g 2 + 2gxg 3 z + g 2 )) + 8c 3 (z 2 - l) (2A/ 2 + g 2 + 2g!g 3 z + g 2 )) 

L( qi ) (96 Cl M 2 (z 2 - 1) ( gl z + g 3 ) + 3c 2 (g!Z + g 3 ) (-20A-/ 2 (z 2 - l) 



61447r 2 F6g 3 (z 2 - l) 2 



+q\ (5 - 2z 2 ) + 6 gl g 3 z + 3g 2 ) - 8c 3 (z 2 - l) (2A/ 2 (3 gi z + g 3 (7 - 4z 2 )) + 3g 3 z 
+g 2 g 3 (4z 2 + 5) + 9g x g 2 z + 3g 3 )) 

L(q 3 ) (96 Cl Af 2 (z 2 - 1) ( qi + g 3 z) + 3c 2 (g! + g 3 z) (-20M 2 (z 2 - l) 



61447r 2 i^gi (z 2 - l) 2 



+3g 2 + 6g!g 3 z + g 2 (5 - 2z 2 )) - 8c 3 (z 2 - l) (2A/ 2 (g x (7 - 4z 2 ) + 3g 3 z) + 3g 3 + 9g 2 g 3 z 

c 2 (g 2 + 2gig 3 z + gf ) 



+g 1 g 2 (4z 2 + 5)+3g 3 z)) 



2048tt 2 F6 (z 2 - 1) 
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n ? ),9 ° A = T^We ~, 2 ^2 J ( 4 : °> "ft' ©J °) c ^ (9? + 2 9193^ + S3) ("4M= (z 2 - l) + q\ + 2 qi q 3 z + q 2 3 ) 
l28F°qiq 3 (z z - 1) 

| c A L{q 2 ) {ql + 2g lg3 z + g 2 ) ((z 2 + 2) (q 2 + 2g l93 z + gg) - 8Af 2 (z 2 - l)) 

30727r 2 F6 g 2 g2 ( z 2 _ !)2 
c 4 £(gi) (-8M 2 (z 2 - 1) ( gl + g 3 z) + g f ( 2 2 + 2 ) + q 2 qaZ ^2 + 5 ) + g^fa 2 + 3ggg) 

30727r 2 F6 gig 2 ( z 2 _ 1) 2 
c 4 L( g3 ) (-8M 2 (z 2 - 1) ( gl z + g 3 ) + 3gfz + 9g 2 iq3 z 2 + g l9 2 z (4z 2 + 5) + gj (z 2 + 2)) 

3072n 2 F%q 2 q 3 (z 2 - l) 2 

c 4 (z (g 2 + ql) + 2 9 ig 3 ) 
30727r 2 F6 gl(?3 (z 2 - 1) ' 



(A.7) 



A Mathematica notebook which contains the above expressions for the structure functions in momentum space is 
available from the authors upon request. 



E. Stephan, S. Kistryn, R. Sworst, A. Biegun, K. Bodek, I. Ciepal, A. Deltuva and E. Epelbaum et al., Phys. Rev. C 82, 
014003 (2010). 

I. Ciepal, S. .Kistryn, E. Stephan, A. Biegun, K. Bodek, A. Deltuva, E. Epelbaum and M. Eslami-Kalantari et al, Phys. 

Rev. C 85, 017001 (2012). 

M. Viviani, L. Girlanda, A. Kievsky, L. E. Marcucci and S. Rosati, arXiv: 1004. 1306 [nucl-th]. 

P. Navratil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007) 
[nucl-th/070i038] . 

D. Gazit, S. Quaglioni and P. Navratil, Phys. Rev. Lett. 103, 102502 (2009) [arXiv:0812.4444l [nucl-th]]. 

R. Roth, J. Langhammer, A. Calci, S. Binder and P. Navratil, Phys. Rev/Lett. 107, 072501 (2011) |arXiv:1105.3T73l 
[nucl-th]]. 

K. Hebeler, J. M. Lattimer, C. J. Pethick and A. Schwenk, Phys. Rev. Lett. 10 5, 161102 (2010) [arXiv:1007.1746l [nucl-th]]. 

G. Hagen, M. Hjorth- Jensen, G. R. Jansen, R. Machleidt and T. Papenbrock. larXiv:12 02.2839 [nucl-th]. 
J. D. Holt, J. Menendez and A. Schwenk, Phys. Rev. Lett. 110, 02 2502 (2013) [arXiv:1207.T509l [nucl-th]]. 
I. Tews, T. Krueger, K. Hebeler and A. Schwenk. la~rXiv:1206.00"25l [nucl-th], 

N. Kalantar-Nay estanaki, E. Epelbaum, J. G. Messchendorp and A. Nogga, Rept. Prog. Phys. 75, 016301 (2012) 
[arXiv:1108. 12271 [nucl-th]]. 

H. -W. Hammer, A. Nogga and A. Schwenk. larXiv: 1210.42731 [nucl-th], 

E. Epelbaum, H. -W. Hammer and U.-G. Meifiner, Rev. M od. Phys. 81, 177 3 (2009) |arXiv:0811. 13381 [nucl-th]]. 
R. Machleidt and D. R. Entem, Phys. Rept. 503, 1 (2011) [arXiv:1105.29T9l [nucl-th]] . 

E. Epelbaum and U.-G. Meifiner, Ann. Rev. Nucl. Part. Sci. 62, 159 (2012) |arXiv:1201.2136l [nucl-th]]. 

D. R. Entem and R. Machleidt, Phys. Rev. C 68, 041001 (2003) [nucl-th/0304018l . 

E. Epelbaum, W. Glockle and U.-G. Meifiner, Nucl. Phys. A 747, 362 (2005) |nucl-th/0405048] . 
U. van Kolck, Phys. Rev. C 49, 2932 (1994). 

E. Epelbaum, A Nogga, W. Glockle, H. Kamada, U.-G. Meifiner and H. Witala, Phys. Rev. C 66, 064001 (2002) 

[nucl-th/0208023] . 

E. Epelbaum, H. Krebs, D. Lee and U.-G. Meifiner, Eur. Phys. J. A 41, 1 25 (2009) |arXiv:0903.T666l [nucl-th]]. 
V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. A 615, 483 (1997) [hep-ph/9611253|. 
N. Fettes, U.-G. Meifiner and S. Steininger, Nucl. Phys. A 640, 199 (1998) |hep-ph /9803266 . 
N. Fettes and U.-G. Meifiner, Nucl. Phys. A 676, 311 (2000) hep-ph/0002162 . 
P. Buettiker and U.-G. Meifiner, Nucl. Phys. A 668, 97 (2000) hep-ph/9908247 . 
J. M. Alarcon, J. M. Camalich and J. A. Oiler, larXiv: 1210.44501 [hep-phj . 
Y.-H. Chen, D.-L. Yao and H. Q. Zheng. larXivTT2 12. 1893 [hep-ph]. 



M. C. M. Rentmeester, R. G. E. Timmermans and J. J. de Swart, Phys. Rev. C 67, 044001 (2003) nucl-th/0302080 



V. Bernard, E. Epelbaum, H. Krebs and U.-G. Meifiner, Phys. Re v. C 77, 06 4004 (2008) arXiv:0712.1967J [nucl-th]]. 

S. Ishikawa and M. R. Robilotta, Phys. Rev. C 76, 014006 (2007) [arXiv:0704.07TT1 [nucl-th]]. 

V. Bernard, E. Epelbaum, H. Krebs and U.-G. Meifiner, Phys. Rev. C 84, 054001 (201 1) [arX"iv:1108.3816l fnucl-thll. 
M. R. Robilotta, Phys. Rev. C 74, 044002 (2006) [Erratum-ibid. C 74, 059902 (2006)] [nucl-th/0610046| . 
J. Golak, D. Rozpedzik, R. Skibinski, K. Topolnicki, H. Witala, W. Glockle, A. Nogga and E. Epelbaum et al., Eur. Phys. 
J. A 43, 241 (2010) [arXiv:0911.4173l [nucl-th]]. 

R. Skibinski, J. Golak, K. Topolnicki, H. Witala, E. Epelbaum, W. Glockle, H. Krebs and A. Nogga et al., Phys. Rev. C 
84, 054005 (2011) [arXiv: 1107. 51631 [nucl-th]]. 

C. Ordonez, L. Ray and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994). 



32 



N. Kaiser, S. Gerstendorfer and W. Weise, Nucl. Phys. A 637, 395 (1998) nucl-th/9802071 . 
H. Krebs, E. Epelbaum and U.-G. Meifiner, Eur. Phys. J. A 32, 127 (20 07) |nucl-th /0703087 . 
N. Kaiser, R. Brockmann and W. Weise, Nucl. Phys. A 625, 758 (1997) [iruc l-th/97060l5p 
R. Machleidt and D. R. Entem, J. Phys. G G 37, 064041 (2010) [arXiv: 1001 09661 [nucl-th]]. 

S. C. Pieper, V. R. Pandharipande, R. B. Wiringa and J. Carlson, Phys. Re v. C 64, 014 001 (2001) [nucl-th/010"2004] . 
H. Krebs, A. Gasparyan and E. Epelbaum, Phys. Rev. C 85, 054006 (2012) I arXiv: 1203 .00671 [nucl-th]^ 
N. Fettes, U.-G. Meifiner, M. Mojzis and S. Steininger, Annals Phys. 283, 273 (2000) [Erratum-ibid. 288, 249 (2001)] 
[hep-ph/000i308 . 



E. Epelbaum, U.-G. Meifiner and J. E. Palomar, Phys. Rev. C 71, 024001 (2005) [nucl-th/04 07037 



W. Glockle, The Quantum Mechanical Few-Body Problem, Springer- Verlag, 1983. 
V. Baru, E. Epelb aum, C. Hanhart, M. Hoferichter, A. E. Kudryavtsev and D. R. Phillips, Eur. Phys. J. A 48, 69 (2012) 

[arXiv: 1202.02081 [nucl-th]] . 

L. Girlanda, A. Kievsky and M. Viviani, Phys. Rev. C 84, 014001 (2 011) [arXiv:1102.4799l [nucl-th]]. 
C. Ordonez, L. Ray and U. van Kolck, Phys. Rev. C 53, 2086 (1996) [hep-ph/9511380|. 
E. Epelbaum, H. Krebs and U.-G. Meifiner, Nucl. Phys. A 806, 65 (2008) |arXiv:0712.1 969 [nucl-th]]. 
E. Epelbaum, H. Krebs and U.-G. Meifiner, Phys. Rev. C 77, 034006 (2008) arXiv:0801.1299 [nucl-th]]. 



